在 解决 实际 问题 的 生产 《或 工程 ) 实践 和 科学 实验 过 程 中 ， 
间 的 函数 关系 来 帮助 我 们 认识 事物 的 内 在 规律 和 本 质 必 性， 
常常 隐 含 在 从 试验 、 观 测 得 到 的 一 组 数据 之 ， 








相对 准确 的 函数 关系 就 成 为 解决 实际 问题 的 关键 。 




















第 十 章 ”数据 拟 合 与 插值 


8 10.1 引言 






































通常 需要 通过 研究 菜 些 变量 之 











而 


这 些 变量 之 间 的 未 知 函数 关系 又 








。 因 此 ， 能 否 根据 一 组 试验 观测 数据 找到 变量 之 间 


例如 在 工程 实践 和 科学 实验 中 ， 常 常 需要 从 一 组 试验 观测 数据 400 了) 门 ， 守 = 0 …) 妖 之 



































中 找到 








来 实现 。 





变量 关 与 攻 
》 一/ (9 的 产生 办 法 











变量 了 之 间 的 函数 关系 ， 一 般 可 用 一 个 近 





以 函数 7 = jz 来 表示 。 函 数 











因 观 测 数据 和 要 求 不 同 而 异 ， 通 常 可 采 











数据 拟 合 与 函数 掺 值 两 种 办 法 




















数据 拟 合 主要 是 考虑 到 观测 数据 受 随机 观测 误差 的 影响 ， 进 而 寻求 整体 误差 最 小 、 能 较 好 


反映 观测 数据 的 近似 函数 多 = /zz) ， 此 时 并 不 要 求 所 





风 =J co) ， 1 = 0, 1 ,7 
函数 插值 则 要 求 近似 函数 》 = Fo) 在 每 一 个 观测 点 X 处 一 定 要 满足 芒 = Joc ) ， 

















1 0 。 在 这 种 | 














诅 到 
可 亚 


































































































中 又 选 

















j 哪 一 种 夺 


















































进行 分 机 ， 尤 为 习 





据 拟 合 虽 然 较 有 效 地 殉 
出 来 的 拟 合 函数 《〈 系 数 ) ， 只 是 拟 合 问题 的 


还 应 该 对 拟 合 函 数 作 区 间 估 计 或 假设 检验 ， 如 果 置 信 区 间 太 大 或 包含 零点 ， 则 由 计算 得 到 的 拟 





























j 拟 合 方法 ， 揪 值 方法 之 中 、 拟 合 方法 












































个 点 估计 ， 还 不 


| 的 近似 函数 了 = 了 (2) 满足 





青 况 下 ， 通 常 要 求 观测 数据 相对 比较 准确 ， 即 不 考虑 观测 误差 的 影响 。 

在 实际 问题 中 ， 通 过 观测 数据 能 否 正确 揭示 某 些 变量 
在 规律 与 本 质 属 性 ， 往 往 取决 于 两 方面 因素 。 其 一 
在 获取 观测 数据 的 过 程 " 








过 间 的 关系 ， 进 而 正确 认识 事物 的 内 
是 观测 数据 的 准确 性 或 准确 程度 ， 这 是 因为 
一 般 存 在 随机 测量 误差 ， 导 致 所 讨论 的 变量 成 为 随机 变量 。 
观测 数据 处 理 方法 的 选择 ， 即 到 底 是 采用 插值 方法 还 是 
直 或 拟 合 技巧 来 处 理 观测 数据 。 折 
问题 则 考虑 了 观测 误差 的 影响 。 但 
情况 下 是 通过 经 验 公 式 获 得 的 ， 



































其 二 是 对 














值 问题 忽略 了 观测 误差 的 影响 ， 而 拟 合 
1 于 观测 数据 客观 上 总 是 存在 观测 误差 ， 而 拟 合 函 数 大 多 数 
因此 要 正确 揭示 事物 的 内 在 规律 ， 往 往 需 要 对 大 量 的 观测 数据 
要 的 是 进行 统计 分 析 。 统 计 分 析 的 方法 有 许多 ， 如 方差 分 析 、 回 归 分 析 等 。 数 
服 了 随机 观测 误差 的 影响 ， 但 从 数理 统计 的 




















和 度 看 ， 根 据 一 个 样本 计算 
能 完全 说 明 其 整体 性 质 。 因 此 ， 





























合 函 数 系数 的 估计 值 就 毫 无 意义 。 这 里 所 采用 的 统计 分 析 方 法 就 是 所 谓 的 回归 分 析 。 另 外 还 可 用 


方差 分 析 的 方法 对 模型 的 误差 作 定量 分 析 。 









































对 于 插值 方法 ， 本 章 在 第 二 节 中 简单 介绍 最 常用 的 插 
于 数据 拟 合 问题 必须 作 区 间 估 计 或 假设 检验 ， 所 以 除了 在 本 章 第 三 节 、 第 四 节 中 介 绎 





























的 数据 拟 合 方法 










































































最 小 二 乘法 的 基本 结论 及 其 Matlab 实 现 问题 外 ， 我 们 在 第 五 节 : 
了 对 数值 拟 合 问题 进行 区 间 估 计 或 假设 检验 的 统计 方法 ， 即 介绍 回归 分 析 方 法 及 其 Matlab 实 现 。 
数据 处 理 问题 通常 情况 下 只 是 某 个 复杂 实际 问题 的 一 个 方面 或 部 分 内 容 ， 
































值 法 的 基本 结论 及 其 Matlab 实 现 问题 。 
最 基本 
专门 介 名 





















































里 所 介 名 





因而 这 

































































的 数据 处 理 方法 一 一 函数 插值 和 数据 拟 合 的 方法 〈 包 括 
部 分 问题 一 一 计算 问题 。 一 般 来 说 ， 对 实际 问题 进行 数学 建 模 需 要 月 























的 情况 下 可 以 单独 使 用 本 章 所 介绍 的 内 容 ， 
建 模 A 题 为 例 说 明 如 何 使 



































用 数值 计算 知识 建立 数学 模型 























故我 们 只 





8$ 10.2 插值 方法 





在 实际 问题 中 所 遇 到 的 插值 问题 一 般 分 为 一 维 插值 问题 和 二 维 



































常用 的 一 维 插值 方法 及 其 一 些 简单 结果 。 























企 本 章 最 



























































也 可 以 是 未 知 的 ) 在 区 间 [w, 世 上 于 二 1 个 互 异 点 屿 处 的 函数 值 了 7 ， 


























sgCo 天 [oO 
维 插值 方法 的 基本 思想 是 ; 
函数 值 蕊 ， 了 = 0 
为 8 的 近似 表达 式 ， 使 得 
JU 六 


然后 计算 AGO 在 区 间 上 [ee 人] 朝 

































































， 求 一 个 足够 光滑 、 





Jj=0， 1，…: 1 。 


交 有 


称 为 插值 区 间 ) 上 点 Y《〈 称 为 撮 

















为 被 插 函 数 ) 在 此 点 的 近似 值 























代数 多 项 式 比 较 简 单 ， 常 用 多 项 式 作为 插值 函 











插值 





多 项 式 的 存在 唯一 


















































数 。 

















简单 便于 计算 的 函 





后 一 节 以 修改 
4， 从 而 解决 实际 问题 的 方法 。 








插值 问题 。 





上 有 若干 阶 导数 ， 如 何 求 出 8C9 在 [@O] 上任 一 点 工 的 近似 值 。 
展 据 g&(x) 在 区 间 [w, 旭 上 于 二 1 个 互 异 点 ”7 ( 称 为 节点 ) 的 
数 /(C50 〔〈 称 为 插值 函数 ) 作 














回归 分 析 ) 通常 上 只 能 解决 实际 问题 中 的 
旧 到 多 方面 知识 ， 上 只 有 很 少 
后 的 美国 91 年 数学 











本 节 主 要 介 


| 








一 维 插值 问题 的 数学 描述 为 ;已 知 某 一 函数 多 = SC0) 〈8(C0 的 解析 表达 式 可 能 十 分 复杂 ， 
= 0.1… 


;”， 还 知道 





《10.1) 

















i 值 点 ) 的 

















。 求 插值 函数 ./ 2) 的 方法 称 为 插值 方法 ，《〈10.1) 称 为 插值 


直 作 为 原 函 数 SCO ( 称 








条 件 。 














假设 了 (Co9) 是 一 个 满足 插值 条 件 〈10.1) 的 次 数 不 超 过 于 的 代数 多 项 式 ， 即 
(10.2) 


xz) =a0 二 QIX 二 十 0 





为 满足 10.1) 的 插值 函数 ， 则 /290 的 4+1 个 待定 系数 0 ， 

















G0 十 和 X0 十 GDX0 十 … 十 


00 十 0 





2 
00 十 01X 十 02X， 十 "… 十 


记 此 方程 组 的 系数 矩阵 为 4 ， 则 








是 范 德 蒙 (Vandermonde) 行 列 式 。 当 707”， 














7 十 个 操 











(10.3 ) 有 唯 解 。. 这 表明 ， 人 外 











2 
2 十 GD2XT 十 十 GT 


大 


下 一 
CjX0 二 .JJ0， 


二 J1， 


基 二 - 
CA 二 水 j， 

















值 节点 Yo ， 人 
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CO 
 “”””) ”满足 











《10.3 ) 





5 光 
“ 互 异 ， 





满足 扣 








惰 互 不 相同 时 ， 此 行列 式 值 不 为 零 。 因 此 方程 组 
ji 值 条 件 〈10.1) 的 

















插值 多 项 式 〈10.2) 存在 唯一 。 
从 几何 上 看 ， 上 述 多 项 式 插值 就 是 过 4+T 个 数据 点 (077)，7 了 -0，1，…， 妖 作 一 条 多 
项 式 曲线 随 三 /90 近似 原 有 曲线 了 = SC) 。 

当 xE [@, 人] 目 < 2 (7 0 用时，8GOD = 7 Co) ， 称 被 插 函 数 8(9D 与 插值 函数 
(多 项 式 ) /0) 之 间 的 差 

RCOD)=SCOOD 一 xD) 

为 插值 函数 〔 多 项 式 ) 7 (0) 的 截断 误差 ， 或 插值 余 项 。 显 然 ， 插 
00 的 次 序 无 关 。 

二 ， 拉 格 朗 日 插值 法 

1 拉 格 衣 日 插值 的 基本 结果 

用 多 项 式 函 数 〈10.2) 作为 插值 函数 时 ， 希 望 通过 解 方程 组 〈10.3) 而 得 到 待定 系数 2o ， 
4 ，…， 色 的 做 法 当 严 比较 大 时 是 不 现实 的 。 方 便 的 方法 是 先 构造 一 组 基 函 数 
1 -=_C 贡 一 DG-3D (一 太 ) 

GCC 0 1 


显然 怀 C0) 是 寻 次 多 项 式 且 满足 


0， 
on 






























































值 多 项 式 的 存在 唯一 与 节点 




















LT 






































































































































了 三 尖 7， 
7 = 六 7 0 





称 五 次 多 项 式 二 (00 为 节点 Xo,2， xzo 上 的 五 次 插值 基 函 数 。 令 


六 = 和 7 

















则 多 项 式 心 C9) 显然 满足 插值 条 件 儿 C) = 六 性 0 7 介 ， 从 而 必 @9) 为 插值 多 项 式 。 插 
值 多 项 式 印 CD) 又 称 为 于 次 拉 格 朗 日 (Lagrange) 插值 多 项 式 ， 册 方程 组 〈10.3) 解 的 存在 叭 
一 性 ，4+ 工 个 互 异 节点 的 于 次 拉 格 朗 日 插值 多 项 式 必 (zx) 存在 唯一 。 

当 8(C0 在 [oO 上 充分 光滑 时 ， 利 用 罗 尔 〈Rolle) 定理 可 推出 ， 对 于 任意 xs [@,0] ， 
插值 多 项 式 妃 Co9) 的 余 项 


RAR COO0)= SO) 一 人 OO) = 















































































































































有 
(+DI! 





Qi (2)， 6 E (0,D) 
OO) 屋 [Tc=->) 
其 ! /0 。 若 可 以 估计 
8 人 EM 
则 可 以 得 到 于 次 拉 格 明 日 插值 的 余 项 估计 


MI ， 
| 
























































实际 上 ， 因 为 妈 周 常 难以 确定 ， 所 以 上 式 并 不 能 给 出 精确 的 误差 估计 。 


















































2 拉 格 衣 日 插值 的 Matlab 实 践 




















Matlab 中 没有 现成 的 拉 格 朗 日 插值 函数 ， 必 须 编写 一 个 M 文 件 实现 拉 格 朗 日 插值 。 



































设 忆 个 节点 数据 以 数组 x0,70 输 入 〈 注 意 Matlab 的 数组 下 标 从 1 开始 ) 
组 兰 输 入 ， 输 出 数组 2》 为 始 个 插值 。 编 写 一 个 名 为 lagrangem 的 M 文 件 : 
function y=1Lagrange (X0vY0vX) ， 
n=1Lendgth (X0) ;m=1Length (X) ， 
for 1=1:m 


























蕊 = 区 ( 工 ) ; 
S=0.0:; 
for K=1:nm 
D=1.0;， 
for ]Jj=1:n 
工 E 了 ~=K 
P=Ppx (Z-X0 (]) )/(x0(K)-x0(]) ) 
end 
ena 
S=DPxyYyO (K)+S7 
enaQ 
Y(I)=S; 
endg 


三 “分 段 线性 插值 法 
1 高 次 插值 多 项 式 的 龙 格 (Runge) 现象 






































用 拉 格 朗 日 插值 多 项 式 包 (29) 作为 区 间 [%, 如 上 连续 函数 SC9 的 近似 函数 ， 在 大 多 数 情况 















































， 了 于 个 插值 点 以 数 





下 ， 忆 (9 的 次 数 越 高 ， 逼 近 8(25) 的 效果 就 越 好 ， 即 误差 | & 429) | 越 小 。 但 是 对 于 高 次 多 项 式 




































































鞭 至 发 生 龙 格 现象 ， 即 在 理论 上 并 不 能 保证 当 寻 趋 于 无 穷 大 时 包 @9) 在 [， 





插值 问题 而 言 ， 往 往 会 造成 插值 多 项 式 如 (25) 的 收敛 性 与 稳定 性 变 差 ， 反 而 逼近 效果 不 理想 ， 











中 上 处 处 收敛 于 














8(x) 。 关 于 这 一 点 ， 在 20 世 纪 初 就 为 龙 格 (Runge) 所 发 现 : 在 三: 上 


用 妈 十 1 个 等 距 节点 作 


被 插值 函数 CO = VGL+ 25x ) 的 插值 多 项 式 卫 CD) ， 则 随 着 地 的 增 大 ， 钱 CD 振荡 越 来 越 大 。 





























ID 











足 条 件 0.726… 到 工 |<1 的 工 ， 忆 (2) 并 不 收敛 于 了 (2 。 












































计算 结果 与 理论 证 明 表明 ， 当 于 趋 于 无 穷 夺 ， 夕 (0 在 区 间 中 部 收 仇 于 (0 ， 但 是 对 于 满 


上 述 例子 说 明 ， 使 用 高 次 多 项 式 插 值 是 危险 的 ， 因 此 在 实际 计算 中 不 能 使 用 高 次 插值 。 这 








反而 启发 我 们 使 用 分 段 插值 方法 ， 即 将 区 间 [2, 人 分 成 一 些小 区 间 ， 在 每 一 个 小 区 间 上 用 低 次 
多 项 式 进行 插值 ， 在 整个 插值 区 间 [,O] 上 就 得 到 一 个 分 段 低 次 多 项 式 插值 函数 区间 的 剖 份 可 
以 是 任意 的 ， 各 小 区 间 上 播 值 多 项 式 的 次 数 的 选取 也 可 按 共 体 问 题 的 要 求 而 选择 。 分 段 低 次 多 
项 式 插值 通常 有 较 好 的 收敛 性 和 稳定 性 ， 算 法 简单 ， 但 插值 函数 光滑 性 变 差 。 第 用 的 分 段 多 项 
式 插 值 法 有 两 类 : 一 类 是 下 面 将 要 介绍 的 分 段 线性 插值 法 ， 另 一 类 是 后 面 将 要 介绍 的 三 次 样 条 















































































































































插值 法 。 
2 分 段 线性 插值 法 
























































假设 区 间 [%, 如 上 的 连续 函数 80) 在 +1 个 节点 4 一 zo < 六 < < 如 二 2 上 的 函数 值 
SC) 二 芒 ，7-0，1，…， 严 。 则 得 到 落 平 面 上 的 靖 +1 个 数据 点 007， 包 ) 。 连 接 相 邻 数据 点 
了) 、(G 得 到 于 条 线段 ， 它 们 组 成 一 条 折线 。 把 区 间 [@,b] 上 这 半 条 折线 段 表示 的 
函数 SG9D 关于 这 慰 十 1 个 数据 点 的 分 段 线性 插值 函数 ， 记 作 7Z(C2) ， 则 7 具有 





















































函数 称 为 被 所 
如 下 性 质 ; 
(DZ7C9 可 以 分 段 表示 ， 在 每 个 小 区 间 [xj-5] 上 ， 它 是 线性 函数 ， 即 
并 一 2 大 一 X 
7 = 了 站 十 2 
IE 贡 一 XiSXS 


(2) 600 三 芒 ， 7 -0， 1，…，7; 
G) 7Z490 在 [2] 上 连续 。 











若 构 造 插值 基 函 数 
一 xze [xu] G= 0 时 全 去 
和 一 Xi 


7 (xz) = ii Xe [xx ] (= 7 时 舍 去 )， 
XX. 

















NA 


革 





0， 其 它 ， 


7(X) 二 》 yo) 
则 2 
人 


而 且 当 &(C9 在 [2,0] 上 一 阶 导数 连续 时 ， 对 于 任意 xs [9, 避 有 


RD HgCO 一 TO 全 放 















































》 


可 2 4 =max{l8 人 


用 7Z40 计算 x 点 的 插值 时 ， 只 用 到 关 左右 的 两 个 节点 ， 计 算 量 与 节点 个 数 姻 二 1 无 关 。 但 刀 
误差 越 小 。 实 际 上 用 数据 点 作 插值 计算 时 ， 分 段 线 性 插值 就 足够 了 ， 如 







































































越 大 ， 分 段 越 多 ， 插 人 
数学 、 物 理 中 用 的 特殊 函数 表 ， 数 理 统 计 中 用 的 概率 分 布 表 等 。 
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3 分 段 线性 插值 的 Matlab 实 现 

用 Matlab 实 现 分 段 线性 插值 不 需要 编制 函数 程序 ，Matlab 中 有 现成 的 一 维 插值 函数 interp1 。 
y=interp1(x0,y0,x,method) 

method 指 定 插值 的 方法 ， 默 认为 线性 插值 。 其 值 可 为 ; 

mearest。 最 近 项 插值 

Tinear ”线性 插值 


'Spline' 立方 样 条 插值 































































































'cubic' 立方 插值 。 
所 有 的 插值 方法 要 求 x0 是 单调 的 。 
当 x0 为 等 距 时 可 以 用 快速 插值 法 ， 使 用 快速 插值 法 的 格式 


为 '#nearestt、'#linear' 、"#Spline' 、'#cubic'。 






































四 “三 次 样 条 插值 法 

1 三 次 样 条 的 基本 结果 

分 段 线性 插值 函数 在 节点 处 的 一 阶 导 数 一 般 不 存在 ， 光 滑 性 不 高 ， 这 影响 了 它 在 诸如 机 械 
加 工 等 领域 〔 希 望 插值 曲线 光滑 ) 中 的 应 用 。 例 如 许多 工程 技术 中 提出 的 计算 问题 对 插值 函数 
的 光滑 性 有 较 高 要 求 ， 如 飞机 的 机 辟 外 形 ， 内 燃 机 的 进 、 排 气门 的 凸轮 曲线 ， 都 要 求 曲线 具有 
较 高 的 光滑 程度 ， 不 仅 要 连续 ， 而 且 要 有 连续 的 曲率 。 在 船舶 .飞机 等 设计 中 ， 面 对 已 知 的 一 些 
数据 点 (oo) ，C 芒 ) ，…，(@ro 力 ) ， 绘 图 员 的 做 法 是 ; 首先 将 这 些 数据 点 描绘 在 平面 
图 纸 上 ， 再 把 一 根 富 有 弹性 的 细 直 条 〈 称 为 样 条 ) 弯曲 ， 使 其 一 边 通过 这 些 数据 点 ， 用 压 铁 因 
定 细 直 条 的 形状 ， 沿 样 条 边沿 绘 出 一 条 光滑 曲线 。 往 往 要 用 几 根 样 条 分 段 完成 上 述 工作 ， 这 时 
应 当 让 样 条 在 连接 处 也 保持 光滑 。 
展 据 力学 理论 进行 分 析 ， 这 样 画 出 的 曲线 ， 在 相 邻 两 数据 点 之 间 实 际 上 是 次 数 不 高 于 3 的 
多 项 式 ， 在 整个 区 间 上 有 连续 的 曲率 。 对 绘图 员 用 样 条 画 出 的 曲线 进行 数学 模拟 ， 导 出 了 样 条 
函数 的 概念 。 

定义 于 区 间 [@,O] 上 的 分 段 函数 C0 ， 若 满足 

1 SC0 在 每 一 小 区 间 芒 六 ? 国 ] 上 是 一 个 三 次 多 项 式 函数 ; 

C) SG9 在 整个 区 间 [@,o] 上 有 连续 的 二 阶 导数 。 
则 称 SGC9 为 [2 加 上 关于 章 分 4 一 xo < 瑟 < < 加 二 2 的 一 个 三 次 样 条 函数 。 

从 而 三 次 样 条 插值 问题 为 ， 给 定 函 数 (40) 在 呈 +1 个 节点 加 ， 国 ，…， 加 得 函数 值 yo， 
力 ，…， 蕊 ， 求 一 个 三 次 样 条 函数 (Cx) ， 使 其 满足 

SG | 

数 汪 00) 称 为 8C9) 的 三 次 样 条 插值 函数 。 
对 于 三 次 样 条 插值 函数 而 言 ， 在 每 一 个 小 区 间 上 有 四 个 待定 参数 ， 共 有 参数 47 个 。 而 在 每 
小 区 间 上 ， 由 插值 条 件 可 得 到 2 个 方程 ， 共 计 27 个 方程 。 光 滑 性 要 求 5(2) 在 每 一 内 部 节点 一 
阶 、 二 阶 导数 连续 ， 从 而 其 左右 导数 相等 ， 这 样 每 一 个 内 部 节点 可 产生 2 个 方程 ， 共 计 “(z 一 卫 
en 为 保证 待定 参数 的 唯 
性 ， 还 差 两 个 方程 。 为 此 ， 常 用 的 方法 是 对 边界 节点 除 函数 值 外 附加 要 求 ， 这 就 是 所 谓 的 边界 
条 件 。 根 据 实 际 问题 的 不 同 ， 三 次 样 条 插值 常用 到 下 列 三 类 边界 条 件 。 
GD m 边 界 条 件 ，2 (@) = 9 (CO) = 了 ， 即 给 定 端点 处 的 一 阶 导数 值 。 由 这 种 边界 条 
件 建立 的 样 条 插值 函数 称 为 SC20) 的 完备 三 次 样 条 插值 函数 。 特别 地 ，》 = 了 "二 0 时 ， 样 条 
线 在 端点 处 呈 水 平 状态 。 
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如 果 Jo 、 功 不 知道 ， 我 们 可 以 要 求 8 (C) 在 端点 处 与 加 、 轧 近 似 相等 ， 即 以 2 交 2， 
为 节点 作 一 个 三 次 拉 格 朗 日 插值 多 项 式 刀 CD ， 以 总 交 --2-3 作 一 个 三 次 拉 格 朗 日 插值 


多 项 式 包 (2) ， 使 得 
9 (a)=N' (al COD)=N(D)， 




































































这 种 边界 条 件 建 立 的 三 次 样 条 插值 函数 称 为 &(z) 的 拉 格 朗 日 三 次 样 条 插值 函数 。 

Ci M 边 界 条 件 ，3 (@) = 0;5 (= 了 ，， 即 给 定 端 点 处 的 二 阶 导数 值 。 特 别 地 
了 = 了 ,二 0 时 ， 称 为 自然 边界 条 件 。 

( 诈 ) 周期 边界 条 件 ， 当 上 = 8(C5) 是 以 2 一 4 二 区 一 20 为 周期 的 周期 函数 时 ， 要 求 SC0) 
也 是 周期 函数 ， 故 端点 处 要 满足 >? (+ 0)=S (0 (+O0)=5 (2 一 0) ， 此 条 件 称 为 周 
期 条 件 。 

给 出 任 一 种 边界 条 件 ， 都 可 以 得 到 两 个 独立 的 方程 这样 42 个 独立 方程 可 唯一 确定 三 次 样 
条 插值 函数 的 47 个 参数 ， 从 而 三 次 样 条 插值 函数 在 理论 上 是 唯一 确定 的 。 

三 次 样 条 插值 函数 的 表达 式 通常 有 两 种 方式 ， 其 具体 的 表现 形式 与 计算 方式 可 查阅 相关 次 
料 。 

当 8(C9) 在 [2] 上 连续 时 ， 满 足 上 述 三 类 边界 条 件 之 一 的 三 次 样 条 插值 函数 *(2) 有 良好 
的 收敛 性 ， 即 
lims(z) = SCO0) 2 全 [a,O] 


























































































































































































































岂 =Imax{x) 一 并 7 】} SO 
其 中 ja 。 而 且 对 于 满足 m 边 界 条 件 或 M 边 界 条 件 的 三 次 样 条 插值 函数 ”2) ， 
当 &(z 在 [如上 四 阶 导 数 连续 时 有 误差 估计 
[网 二 8 亲属 BE 有， 
3 1 3 


co = 一 cl= 一 c = 二 二 人 
人 
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2 三 次 样 条 插值 的 Matlab 实 现 

在 Matlab 中 数据 点 称 之 为 断 点 。 如 果 三 次 样 条 插值 没有 边界 条 件 ， 最 常用 的 方法 ， 就 是 条 
j 非 扭 结 〈notra-knot) 条 件 。 这 个 条 件 强迫 第 1 个 和 第 2 个 三 次 多 项 式 的 三 阶 导数 相等 。 对 最 后 
一 个 和 倒数 第 2 个 三 次 多 项 式 也 做 同样 地 处 理 。 

Matlab 中 三 次 样 条 插值 也 有 现成 的 函数 : 

y=interp1(x0,y0,x,Spline'); 




















































































































y=Spline(X0,y0,X); 

pp=csape(X0,y0,conds)， 
pp=csape(X0,y0,conds,valconds)，y=ppval(pp,x)。 
x0,y0 是 已 知 数据 点 ，x 是 插值 点 ，y 古 插值 点 的 函数 值 。 
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对 于 三 次 样 


和 


数 


亡 wv 

















"complete' 
则 按 缺 省 情况 处 
mot-a-knot 








perliodic' 


'Second' 





条 插值 ， 我 们 提倡 使 
直 ， 必 须 调用 函数 ppval。 
pp=csape(x0,y0): 使 用 默认 的 边界 条 伯 
pp=csape(x0,y0,conds,valconds) 中 的 conds 指 定 扣 
边界 为 一 阶 导数 ， 


理 。 



































如 














函数 csape，csape 的 返回 














值 





1 














非 扭 结 条 件 


周期 条 件 





二 阶 导 数 的 缺 省 值 为 [0, 0]。 
variational， 设置 边界 的 二 阶 导 数值 为 [0.0]。 








个 元 素 表示 右边 界 的 条 作 








valconds 给 出 。 


























详细 情况 请 
例 1 机 床 加 


























使 
























































边界 为 二 阶 导 数 ， 二 阶 导数 的 值 






































用 帮助 help csape。 
























































值 是 pp 形式 ， 要 求 插值 点 的 函 

















F， 即 Lagrange 边 界 条 件 。 
的 边界 条 件 ， 其 值 可 为 : 
阶 导数 的 值 在 valconds 参 数 中 给 出 ， 若 忽略 valconds 人 参数 ， 




















在 valconds 参 数 中 给 出 ， 若 忽略 valconds 参 数 ， 





对 于 一 些 特殊 的 边界 条 件 ， 可 以 通过 conds 的 一 个 1X2 矩阵 来 表示 ，conds 元 素 的 取 值 为 0， 





conds(D)=j 的 含义 是 给 定 端点 的 了 阶 导数 ， 即 conds 的 第 一 个 元 素 表示 左边 界 的 条 件 ， 第 二 











F ， conds=[2,1] 表 示 左 边界 是 二 阶 导 数 ， 二 


一 组 数据 4z, 攻 ) 给 出 〈 在 平面 情况 下 ) ， 用 程 



































边界 是 一 阶 导 数 ， 对 应 的 值 






























































待 加 工 零件 的 外 形 根据 工艺 要 求 | 控 铣 床 加 
工时 每 一 刀 只 能 沿 并 方向 和 了 方向 走 非常 小 的 一 步 ， 这 就 需要 从 已 知 数据 得 到 加 工 所 要 求 的 步 
长 很 小 的 42y) 坐标 。 

表 中 给 出 的 风 了 数据 位 于 机 辟 断 面 的 下 轮廓 线 上 ， 假 设 需要 得 到 站 坐标 每 改变 0.1 时 的 了 
坐标 。 试 完成 加 工 所 需 数据 ， 画 出 曲线 ， 并 求 出 x= 0 处 的 曲线 斜率 和 13<xS15 范围 内 上 的 
最 小 值 。 

YX 035791l1213 14 15 
) 01412172.02.12.01.812 1.0 1.6 
要 求 用 Lagrange、 分 段 线 性 和 三 次 样 条 三 种 插值 方法 计算 。 
解 ”编写 以 下 程序 ; 
XO= [0 3 D 7 9 工 2 下 14 15]，; 
溺 0 二 1 于 5 人 下 0 0 8 | 
X=0:0.1:15)， 
yl=1agrange (x0,y0v,x); s 前 面 编 写 的 拉 格 朗 日 插值 函数 





yY2=1interp1(X0v,y0yX) ， 














y3= 衬 Ti 七 全 荆 加 工 (区 0 707 广 7 “SILnEG7 ) ， 


PP1=csape(X0vyYy0) ; 


Y4=pPPVvVal (PP1，X) ， 


PP2=csape (X0,Yy0， seconad' ) ， 





Y5=PPVal (PP2，X) ; 
[x'vyl'vy2'vy3'vy4'vyY5 ] 
Subp1lot (27271) 

节 生 9 芋 (其 0 庆 训 0 下 元 区 六 立 芋 ) 
titlIe(' Iagzange  ) 
Subp1lot (27272) 

这 王 @ 丰 (区 D 关 学 97 下 7 部 公 ) 


title('Piecewise Linear') 








Subp1lot (27273) 

全 工 G 蕊 (区 07 0 和 和 生 3) 
七 斌 七 下 全 和 下 业 革 让 训 工 ) 
Subp1lot (27274) 

生 1O 蕊 (区 0 7 训 0 二 77 区 7 立 4) 
七 i 廿 Je ('SPLine2 ') 
QX=Qiff (xX) ， 
Qy=dqiff(y3) ， 
Qy_ _qx=qy./dx'， 
Qqy qxu=dqy Qqx (1) 
ytemp=y3 (131:151) ， 





ymin=min (ytemp) ， 
indqex=finda(y3==ymin) ， 
Xmlin=X(inadqex) ， 
[xminyyminl] 


计算 结果 略 。 








可 以 看 出 ， 拉 格 衣 日 插值 的 结果 根本 不 能 应 用 ， 分 段 线 性 插值 的 光滑 性 较 差 〈 特 别 是 在 
X=]T4 附近 弯曲 处 ) ， 建 议 选 用 三 次 样 条 插值 的 结果 。 










































































五 “一 维 插值 总 结 





插值 函数 一 般 是 已 知 函数 的 线性 组 合 或 者 称 为 加 权 平 均 。 在 已 知 数据 点 较 少 时 ， 插 值 技术 
在 工程 实践 和 科学 实验 中 有 着 广泛 而 又 十 分 重要 的 应 用 。 例 如 在 信息 技术 中 的 图 像 重 建 、 图 像 放 
大 过 程 中 为 避免 图 像 失 真 、 扭 曲 而 增加 的 插值 补 点 ， 建 筑 工 程 的 外 观 设 计 ， 化 学 工程 试验 数据 











与 模型 分 析 ， 天 文 观测 数据 、 地 到 
术 的 应 用 是 不 可 或 缺 的 。 

































































信息 数据 的 处 理 ， 社 会 经 济 现象 的 统计 分 析 等 方面 ， 插 值 技 
































插值 技术 《或 方法 ) 远 不 止 这 里 所 介绍 的 这 些 ， 但 在 解决 实际 问题 时 ， 对 于 一 位 捐 
一 














值 问题 











而 言 ， 前 面 介绍 的 插值 方法 已 经 足够 了 。 剩 下 的 问题 关键 在 于 什么 情况 下 使 用 、 怎样 使 用 和 使 用 

















何 种 插值 方法 的 选择 上 。 





拉 格 衣 日 插值 函数 在 整个 插值 区 间 上 有 统一 的 解析 表达 式 ， 其 形式 关于 节点 对 称 ， 光 滑 性 
好 。 但 缺点 同样 明显 ， 这 主要 体现 在 高 次 插值 收敛 性 差 〈《 龙 格 现象 ) ;增加 贡 点 时 前 期 计算 作 
废 ， 导 致 计算 量 大 ;一 个 节点 函数 值 的 微小 变化 〈 观 测 误差 存在 ) 将 导致 整个 区 间 上 插值 函数 









































都 发 生 改变 ， 因 而 稳定 性 差 等 几 个 方面 。 因 此 拉 格 朗 日 插值 法 多 用 于 至 
日 插值 方法 进行 插值 计算 时 通常 选取 有 < 7 。 







































































8 论 分 析 ， 在 采 j 





日 拉 格 衣 














分 段 线性 插值 函数 〈 仅 连续 ) 与 三 次 样 条 插值 函数 〈 二 阶 导数 连续 ) 虽然 光滑 性 差 ， 但 他 

















们 都 到 服 了 拉 格 朗 日 插值 函数 的 缺点 ， 不 仅 收敛 性 、 稳 定性 强 ， 而 且 方 法 简单 实 月 








因而 应 用 十 分 广泛 。 
六 ”二 维 插 值 简介 



































日 ， 计 算 量 小 。 


前 面 讲 述 的 都 是 一 维 插值 ， 其 特点 是 节点 为 一 维 实数 组 ， 插 值 函数 是 一 元 函数 〈 曲 线 ) 。 
若 节点 是 二 维 数组 ， 则 插值 函数 就 是 二 元 函数 ， 即 曲面 。 例 如 在 某 区 域 测量 了 若干 点 《节点 ) 




































































的 高 程 〈 节 点 值 ) ， 为 了 画 出 该 区 域 较 精确 的 等 高 线 图 ， 就 要 先 插 入 更 多 的 点 《插值 点 ) 























算 这 些 点 的 高 程 〈 插 值 ) ， 这 就 是 一 个 二 维 插值 问题 。 






























































;, 强 


二 维 插值 问题 的 数学 描述 为 ;已 知 二 元 函数 go 尹 在 某 矩 形 区 域 尽 =[@,O]x[cd] 上 所 
异 节点 5 了/ 的 函数 值 宝 ， 如 何 求 出 在 及 上 任 一 点 Co 切 处 的 函数 值 &Cxc,y) 的 近似 值 。 






























































一 维 插值 方法 的 基本 思想 是 ， 根 据 已 知 数据 点 (25727) ， 构 造 一 个 具有 一 定 光 滑 性 、 简 














单 易于 计算 的 函数 2= /02J) 〔 已 知 曲面 ) 作为 = 8(C2y) 〔( 未 知 曲 
?= Ce 力 通 过 已 知 数据 点 ， 即 了 Go) 一列 


SG 切 的 近似 值 。 












































掉 ) 的 近似 ， 使 曲 
































常见 的 二 维 插值 可 分 为 网 格 节 点 捐 


























































































































田 








。 然 后 计算 点 Co 切 的 函数 值 / Co 攻 ) 作为 


直 和 散乱 节点 插值 。 其 中 网 格 节点 插值 适用 于 节点 比较 


规范 的 情况 ， 即 在 包含 所 给 节点 的 矩形 区 域内 ， 节 点 由 两 族 平 行 于 坐标 轴 的 直线 的 交点 所 组 成 。 




















散乱 节点 插值 适用 于 一 般 的 节点 ， 多 用 于 节点 不 太 规 范 《〈 即 节点 为 两 族 平行 于 坐标 轴 的 直 























部 分 交点 ) 的 情况 。 
1 网 格 节点 插值 法 
网 格 节点 插值 问题 可 简 述 为 ， 已 知 数据 点 (27) ， 其 中 写 -1 





























任 全 2) 访 处 的 插值 z 。 
网 格 节点 插值 法 主要 有 以 下 几 种 形式 : 


(1) 分 片 线性 插值 





























1 Q= 为 <X< <X = CC= 力 <J< < =C 求 R 上 人 














ee 171 ， 7 -1 























[ 线 的 


，2， 


F 插值 点 (xz,J) 


分 片 线性 插值 即 是 在 以 25， 切 ) (Co ， (Goal JP) Go 7 为 顶点 的 小 矩形 他 





上 作 线 性 插值 。 其 分 片 插值 函数 为 : 




















JJj4+l 一 


J 和 了 革 人 X) 十 也 
当 Vi 大 冬 Ki+H1 。 Xi+Hl 一 时 ， 
JoJ 攻 =21+(2H 一 2 认 六) 十 (Za 一 ZJ 一) 
了 jH 少 / 


| 
J 0 切 = 巧 二 (2 一 2 一 六) 十 (25 一 37) 一 攻 ) 
显然 ， 分 片 线性 插值 函数 在 尽 上 连续 。 
CO) 分 片 双 线性 插值 
分 片 双 线性 插值 曲面 由 一 片 一 片 空间 二 次 曲面 构成 ， 其 分 片 函数 表达 式 为 ; 
Joom)=(Caz+BJCc+D)， Cos RN 
中 四 个 待定 系数 可 由 如 四 个 顶点 值 唯一 确定 。 但 分 片 双 线性 插值 曲面 在 及 上 不 连续 。 
G) 分 片 双 三 次 样 条 插值 
分 片 双 三 次 样 条 插值 函数 在 每 一 片 〈 即 小 矩形 &o ) 上 的 表达 式 为 
71)=(4+4Xx+40 +4)(B+PDXT+TBX 二 4 ) 
其 中 待定 系数 省 ， 忌 人 = 2.3,4) 可 由 各 四 个 顶点 的 值 及 插值 了 攻 Ge 功 在 和 了 方向 的 
滑 性 “ 即 庆 、 廊 、 户 、 和 连续 ) 和 相应 的 边界 条 件 所 唯一 确定 。 
2 散乱 节点 插 值 
散乱 节点 插值 问题 可 简 述 为 ， 已 知 函数 三/ (2 切 在 矩形 域 尽 上 散乱 分 布 的 个 互 异 节 
点 (5* 世 7 处 的 函数 值 芍 ， 求 及 上 任 一 插值 点 (2J) 性 0 六 的 插值 z 。 
常用 的 方法 是 反 距 离 加 权 平 均 法 ， 或 称 Shepard 方 法 。 其 基本 思想 是 ， 在 点 CD 
公 C5 了 ， 定 义 其 插值 函数 的 函数 售 为 节点 然 函数 值 按 C5 了 ) 与 节点 距离 的 某 种 形式 反比 
作为 权重 的 加 权 平均 。 例 如 车 记 万 二 VCx 一 芋 ) 十 作 一 蕊 )》 ， 则 插值 函数 〈 曲 面 ) 可 定义 为 
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庆 方 = 0， 


Jo DJ)= (2J)27 万 夫 0， 


殉 ( 芒 = 记 并 二 


尼 / 分 交 

这 样 定义 的 插值 曲面 是 全 局 相关 的 ， 对 曲面 上 任 一 点 作 数 值 计 算 都 要 涉及 到 全 体 已 知 数据 ， 
因而 在 已 知 数据 量 大 的 情况 下 计算 量 相当 大 ， 而 且 插值 曲面 在 节点 (25? 世 ) 处 不 光滑 。 虽 然 如 
此 ， 由 于 该 方法 思想 简单 ， 故 有 种 种 改进 方法 。 














































































































3 二 维 插值 的 Matlab 实 现 
Matlab 中 有 一 些 计算 二 维 插值 的 程序 。 如 
Z=interp2(X0,y0,z0,x,y,method) 




















其 
组 ， 


















































x0,y0 分 别 为 丈 维和 于 维 向 量 ， 表 示 节 点 ，z0 为 允 X 丸 维和 矩阵 ， 表 示 节 点 值 ，x,y 为 一 维 数 
表示 插值 点 ，x 与 y 应 是 方向 不 同 的 向 量 ， 即 一 个 是 行 向 量 ， 另 一 个 是 列 向 量 ，z 为 矩阵 ， 



































表示 得 到 的 插值 ，method' 的 用 法 同上 面 的 一 维 插值 。 





其 


























如 果 是 三 次 样 条 插值 ， 可 以 使 用 命令 
pp=csape({x0,y0},z0,conds,valconds)，y=ppval(pp,{xy 放 

x0,y0 分 别 为 到 维和 于 维 向 量 ，z0 为 亚 关 于 维 矩阵 ， 有 具体 使 用 方法 同一 维 插值 。 

例 2 在 一 丘陵 地 带 测量 高 程 ， 守 和 ?方向 每 隔 100 米 测 一 个 点 ， 得 高 程 如 下 表 ， 试 拟 合 一 











































































































曲面 ， 确 定 合 适 的 模型 ， 并 由 此 找 出 最 高 点 和 该 点 的 高 程 。 


























次 


了 100 200 300 400 500 
100 636 697 624 478 450 
200 698 712 630 478 420 
300 680 674 598 412 ”400 
400 662 626 552 334 310 


解 ”编写 程序 如 下 : 














ClLearyclc 














X=100:100:5007 

y=100:100:400， 

Z=[636 697 624 478 450 
698 712 630 478 420 
680 674 598 412 400 
662 626 避 号 儿 334 3101]， 


PP=csape({xv yjvz ") 


人 有 


=fnval (Pb {(100:10:500，100:10:4001】) 








CZ2=1nterp2 (xyryrzy100:10:500，[100:10:400]7" spPline') 


8 3 数据 拟 合 
在 科学 计算 中 经 常 要 建立 实验 数据 的 数学 模型 。 给 定 函 数 的 实验 数据 ， 需 要 用 比较 简单 和 





合适 的 函数 来 逼近 〔〈 或 拟 合 ) 实验 数据 。 这 种 逼近 的 特点 是 : 





Ca) 适度 的 精度 是 需要 的 ; 
Cb) 实验 数据 有 小 的 误差 

(@ 对 于 某 些 问 题 ， 可 能 有 某 些 特殊 的 信息 能 够 用 来 选择 实验 数据 的 数学 模型 。 

逼近 离散 数据 的 基本 方法 就 是 曲线 拟 合 ， 常 采用 最 小 二 乘 拟 合 。 

曲线 拟 合 问题 的 数学 描述 是 ， 已 知 一 组 〈 二 维 ) 数据 45， 芒 ) ，7 一 12 于 《〈 即 平面 上 的 

































































天 个 点 C 态 ) = 2 天) ， 瑟 王 不 相同 。 寻 求 一 个 函数 〈 曲 线 ) 了 = /Co) ， 使 7 在 
某 种 准则 下 与 所 有 数据 点 最 为 接近 ， 即 曲线 拟 合 得 最 好 。 

















一 “线性 最 小 二 乘 拟 合 
线性 最 小 二 乘法 是 解 次 曲线 拟 合 问题 最 常用 的 方法 ， 基 本 思路 是 ， 令 
JG0=aiGCOD)TG2D(OOD 二 十 Go (10.4) 
其 中 族 COj= Cn< 六 是 一 组 事先 选 定 的 线性 无 关 的 函数 ，{oxeje-: 是 一 组 待定 系数 。 寻 求 系 
数 1oeji 使 得 六 与 G5) 的 距离 性 =12,…) 的 平方 和 最 小 。 这 种 准则 称 为 最 小 二 乘 准 
则 ， 其 求 系数 {oej = 的 方法 称 为 线性 最 小 二 乘 拟 合 方法 。 
1 系数 tej=: 的 求法 










































































若 记 
Jan)=》 妇 =》 [co)-7 
已 1 i=1 ， 《10.5 ) 
则 了 为 和 4 的 二 次 函数 。 由 数学 分 析 〔〈 或 高 等 数学 ) 的 极 值 理 论 ，v 达到 最 小 的 充分 必要 



































97 _ 


条 件 是 0…Qawn 满足 9 (三 1 zz2) 。 于 是 得 到 求 使 了 达到 最 小 的 和 ，.…，9w 的 方法 
是 求解 线性 方程 组 〈 称 为 法 方程 组 ) 


和 PCo)I onco)-o=0 (7 =1 ,70) 


》 




















即 求解 线性 方程 组 
ax[2 六 GO]= 记 六 OCT=L 
大 =1 大 =1 冲 





《10.6) 


10) 7 
R=| : : 
站 (Co) Cj 4=[ooeao] 了 =(O0 





则 方程 组 〈10.6) 可 表示 为 
R R4= 及 了。 (10.7) 


由 于 当世 (0 (让 线性 无 关 时 ， 玉 列 满 秩 ， 及 7 及 可 道 ， 所 以 方程 组 〔〈10.7) 有 唯一 














4=(RR)RT7 





用 《〈10.4) 作 线 性 最 小 二 乘 拟 合 的 误差 通常 考虑 以 下 两 种 形式 ; 
WE 

了 最 小 平方 误差 : 生 ; 

最 大 偏差 ， 0， = max | J 9 


1<i<7 














2 函数 组 巡 CODj=; 的 选取 
面 对 一 组 数据 (2 六 因 三 上 2 于， 用 线性 最 小 二 乘法 作曲 线 拟 合 时 ， 首 要 的 、 也 是 关键 






































的 一 步 是 恰当 地 选取 五 (5) 六 (5) 。 如 果 通 过 机 理 分 析 ， 能 够 知道 了 与 并 之 间 应 该 有 什么 样 























的 函数 关系 ， 则 五 (7 zw CD) 容易 确定 。 若 无 法 知道 与 并 2 
Go 六) 12… 浊 作 图 ， 直 观 地 判断 应 该 用 什么 样 的 














(i) 直线 了 = 二 0X 十 22 ; 

















间 的 关系 ， 通 常 可 以 将 数据 
































线 去 作 拟 合 。 人 们 常用 的 曲线 有 





























(ii) 多 项 式 了 = 0X 二 十 Gonx+Qor (一般 本 =2.3， 不 宣 太 高 ) ， 














了 0 = 1] 
多 、 2 se 扎 ! 、 QQ AAA 
《证 ) 双 曲 线 《〈 一 支 ) 和 ， 拟 合 前 作 变 量 蔡 换 工 求解 “2 较 简单 ， 
































2 GD 区 1 三 
Civ) 指数 曲线 “了 = aie” ， 拟 合 前 作 变 量 代 换 z= 世 》，“ 一 /zx ， 则 指数 曙 线 
》 二 4 ”转化 为 关于 下 4， 42 的 线性 函数 了 二 下 4 + 42f ， 这 样 做 求解 4142 较 简单 

















在 实际 计算 过 程 中 ， 面 对 一 组 已 知 数据 ， 至 








j 底 用 什么 样 的 


























线 拟 合 最 好 ， 可 以 在 直观 判断 



































的 基础 上 ， 选 几 种 曲线 分 别 拟 合 ， 然 后 比较 ， 看 哪 条 曲线 的 最 








二 ， 非 线性 最 小 二 乘 拟 合 
非 线性 最 小 二 乘法 是 假设 /7 (0) 是 待定 系数 
求 其 系数 {oejk- 。 








fc 和 的 任意 








小 二 乘 指标 最 小 。 


非 线 性 函数 ， 在 最 小 二 乘 准则 下 








例如 上 述 人 们 和 常 用 的 双 曲 线 和 指数 曲线 就 是 非 线性 最 小 二 乘 拟 合 中 最 常用 的 非 线性 函数 ， 











只 不 过 在 上 面 使 用 中 将 它们 转变 成 线性 最 小 二 乘 拟 合 方法 。 











对 于 给 定 的 实验 数据 ， 通 常 应 根据 实验 数据 的 走向 、 趋 势 选 择 合适 的 数学 模型 ， 即 拟 合 函 

















数 。 例 如 当 实 验 数 据 具 有 单调 性 和 凸 性 时 ， 可 选择 下 述 适当 














的 数学 模型 了 = 了 (9) 来 拟 合 实验 数 


JCD=ae”， JCOD=ae JCD=axr xz)=C+D/X 











等 ， 其 中 4 、2 为 参数 。 








在 有 可 能 的 情况 下 ， 一 般 将 非 线性 拟 合 函数 转化 轴线 性 朱 合 函数 求解 ， 这 一 方面 是 如 此 求 
解 简单 ， 另 -方面 也 是 因为 一 般 情况 下 求解 法 方程 组 9 
(&592，…Qx) 通常 仅 是 了 的 驻 点 ， 不 一 定 是 极 值 点 。 也 可 以 直接 解 了 极 小 化 问题 。 



































三 “最 小 二 乘 拟 合法 的 Matlab 实 现 
1 解 方程 组 方法 


在 上 面 的 记号 下 ， 
Joan)=|R4-Z| 











Matlab 中 的 线性 最 小 二 乘 的 标准 型 为 
Min |R4- 习 








和 
2 ， 


命令 为 4=R\7 。 


PT 




















多 = 四 得 到 的 


























本 2 
例 3 用 最 小 二 乘法 求 一 个 形 如 了 二 4 二 Ox 的 经 验 公 式 ， 使 它 与 下 表 所 示 的 数据 拟 合 。 











xx 19 253 3l1 38 44 


了 


19.0 32.3 49.0 73.3 97.8 


解 ”编写 程序 如 下 


X= [19 
y=[19.0 





2 


25 


3 


31 38 44] "， 
49.0 风 33 371.8j 


[Onegs (57 IT) 7。 2] 7 


ab=LzANY 


XO=19:0.1:44; 


yY0=ab (1)+ab(2)r*x0.^ 人 2; 


薄 王 G 巧 ( 交 7 训 六 商 7 瘙 OU 芋 昌 
2 线性 最 小 二 乘 拟 合 《多项式 拟 合 ) 方法 


在 线性 最 小 二 乘 拟 合 中 ， 月 
人 (CC)} 一 人 克 生 汪 交 








YI 
< 





? 


a=polyfit(x0,y0,m)， 
其 中 输入 参数 x0,y0 为 要 拟 合 的 数据 ，m 为 拟 合 多 项 式 的 次 数 ， 输 出 参数 a 为 拟 合 多 项 式 
y=amxn+…+alx+a0 系 数 a=[ aa，…，ali，ao]。 
多 项 式 在 x 处 的 值 y 可 用 下 面 的 函数 计算 
y=polyval (a, X) 。 





















































例 4 某 乡镇 企业 1990-1996 年 的 生产 利润 如 下 表 : 


年 份 
利润 〈 万 元 ) 


试 预 调 1997 刍 











1990 1991 1992 1993 1994 
70 122 144 152 174 
FE 和 1998 年 的 利润 。 


解 ” 作 已 知 数据 的 的 散 点 图 ， 


X0O= [1990 


yY0=[70 


1991 


122 


1 。 
于 


的 较 多 的 是 多 项 式 拟 合 。 如 果 取 
日 始 次 多 项 式 拟 合 给 定数 据 ， 则 Matlab 中 有 现成 的 函数 








1995 1996 
196 202 


二 992 工 993 -上 994， 荆 995，- 工 996] 


144 152 174 196 


ELo 蕊 (区 0 YY0 5xT) 


发 现 该 乡镇 企 、 
数 来 预测 该 乡镇 企 ， 


X0O= [1990 


Y0O= [70 





上 的 年 生产 利润 几乎 直线 上 升 。 






































工 991 


122 


202] 


因此 ， 我 们 可 以 用 了 三 ix 20 作为 拟 合 函 








上 未 来 的 年 利润 。 编 写 程序 如 下 : 


1992 1993 1994 1995 1996]:， 


144 22 174 196 


a=DPolLyYyfit(x0vy0y1) 
(ay1997) 


Y97=pPo] 


YVa] 








Y98=Ppo ] 


YVa] 





(av 





二 998)) 





2021]， 


求 得 和 二 20，9o 二 0705X10 ，1997 年 的 生产 利润 997=233. 4286，1998 年 的 生产 利 





润 y98 


=253.9286。 


3 非 线性 最 小 二 乘 拟 合 
Matlab 的 优化 工具 箱 中 提供 了 两 个 求 非 线性 最 小 二 乘 拟 合 的 函数 : curvefit 和 1leastsq。 使 




















] 


CUrvefit 












































这 两 个 命令 时 ， 都 要 先 建立 M 文 件 fun m， 但 它们 定义 了 (5) 的 方式 是 不 同 的 。 





设 已 知 xdata=(xdatal, xdataz,…, xdata2 ),ydata=(ydatal, ydataz ,ydata7 )， 


CUrveTf it 









































以 求 含 参量 x〈 向 量 ) 的 向 量 值 函 数 
F(x, xdata)=(f(x, datal)，…,f(x, xdata2 ))TI 




















中 的 参 变 量 x《〈 向 量 ) ， 使 得 


最 小 。 








也 


i=1 ( 


输入 格式 为 : 














| 
| 


(1) x=curve1 
《2) x=curve1 
(3) x=curvel 
(4) [x，opti 
(5) [x，opti 
(6) [x，opti 
输出 目标 函数 值 格 式 : 作 fan(x,xdata). 
x0 为 从 代 初 值 ，options 为 控制 参数 。 

2 
设 已 知 





leastsq 














leastsq 用 以 求 含 


中 的 参 


最 小 。 








F(x, xdatai)-ydatai)“ 


fit (fun', x0, xdata, ydata) ; 

Fit (fun', x0, xdata, ydata, options) ; 

Fit (fun', x0, xdata, ydata, options，'erad) ; 
onsj=curvefit (fun', x0, xdata, ydata, …) ; 


ons, funvalj=curvefit (fun', x0, xdata, ydata, …) ; 





ons, funval, Jacob]=curvefit (fun', x0, xdata, ydata,…) . 




















xdata=(xdatal, xdataz …, xdata"),ydata=(ydatal, ydataz ,ydata7 )， 























参量 x〈 向 量 ) 的 向 量 值 巩 数 
fx)=(fGx)，f2(x) ff (Cx)) 
变量 x〈 回 量 ) ， 使 得 


























f(x) 节 (x)= 和 (2 f(x) + + (COx) 


其 中 


fi (x)=F(x, xdataiy ydatai) 
输入 格式 为 : 
(1) x= leastsq (fun', x0) ; 


(2) x= leastsq (fun', x0, options) ; 


(3) x= leastsq (fun', x0, options，'egrad) ; 


(4) [x，options]= leastsdqd (fun', x0,…) ; 


(5) [x，options, funval]j= leastsq (fun', x0,，…) ; 
-0.02 乳 


例 5 用 下 面 一 组 数据 拟 合 函数 c() = 4+ pe ”中 的 参数 光 六 A 。 




















人 100 200 300 400 300 000 700 800 900 1000 


, X103; 
7 0 4.34 4.99 $.3$ .63 .90 0.10 60.26 6.39 6.$0 0.59 





解 ”该 问题 即 阶 最 优化 门 是 
min 严 (a,D, 丰 ) = 志 二 De 一 c 中 
1 用 命令 curvefit。 此 时 
RE 
G) 编写 M 文 件 curvefunl.m 
function 仁 curvefun1(x,tdata) 
传 x(])+x(2)#exp(-0.02*x(3)*tdata) % 其 中 x()=a x(2)=b; x(3)=Kk; 
(2) 输入 命令 
tdata=100:100:1000 
cdata=le-03*#[4.54,4.99,$.33,$.63,5.90,6.10,6.26,6.39,6.30,6.59]; 


x0=[0.2.0.05,0.005]; 





X=CuUrvefit( curvefun1l 7,x0,tdata,cdata) 
仁 curvefun1(x,tdatal) 
《3) 运算 结果 为 : 
Xx=0.0070 -0.0030 0.1012 
仁 





Columns 1l through 7 
0.0045 0.0050 0.0054 0.0057 0.0059 0.0061 0.0063 
Columns 8 through 10 
0.0064 0.0065 0.0066 
即 拟 合 得 a=0.0070，b=-0.0030，k=0.0066。 

2 用 命令 leastsq。 此 时 

人 
Xx=(a3，b，Jk) 

G) 编写 M 文 件 curvefun2.m 

function 仁 curvefun2(X) 


tdata=100:100:1000; 


cdata=le-03*[4.34,4.99,3.33,3.603,3.90,6.10,0.26,6.39,6.30,0.59]; 
仁 cdata-x(]1)-x(2)*exp(-0.02*x(3)*#tdata) % 其 中 x(1)=a; x(2)=b:; x(3)=k; 
(2) 输入 命令 


x0=[0.2.0.05,0.005]: 














X=leastsq( curvefun2”,x0) 
仁 curvefun2(X) 
(3) 运算 结果 为 : 
Xx=0.0070 -0.0030 0.1012 
仁 1.0e-005S* 
Columns l through 7 
0.0221 0.2081 -0.3933 -0.2872 0.2973 0.3561 0.0693 





Columns 8 through 10 
-0.2327 -0.0970 0.0296 
可 以 看 出 ， 两 个 命令 的 计算 结果 是 相同 的 。 














84 函数 的 最 小 二 飞盘 近 

在 科学 计算 中 我 们 还 常 遇 到 已 知 函 数 的 逼近 问题 。 由 于 电子 计算 机 只 能 做 算术 运算 ， 因 而 
在 计算 机 上 计算 数学 函数 〈 如 ”等 在 有 限 区 间 上 的 计算 ) 必须 用 其 它 简 单 的 函数 如 多 项 式 或 
有 理 分 式 ) 来 逼近 ， 且 用 该 函数 代 蔡 原来 精确 的 数学 函数 进行 计算 。 这 种 函数 逼近 的 特点 是 : 

Ca) 要 求 是 高 精度 吉 近 ， 

Cb) “要 快速 计算 〈 计 算 量 越 小 越 好 ) 。 

前 面 讲 的 曲线 拟 合 是 已 知 一 组 离散 数据 {2 世 )7 三 2 录 ， 选 择 一 个 较 简单 的 函数 
VC0 ， 如 多 项 式 ， 在 一 定 准则 如 最 小 二 乘 准则 下 ， 最 接近 这 些 数据 。 

如 果 已 知 一 个 较为 复杂 的 连续 函数 了 zx) xs [2 ， 要 求 选择 一 个 较 简单 的 函数 7 (2 ， 








































































































































































































在 一 定 准则 下 最 接近 (2) ， 就 是 所 谓 函 数 逼 近 。 
与 曲线 拟 合 的 最 小 二 乘 准 则 相对 应 ， 函 数 逼 近 常 用 的 一 种 准则 是 最 小 平方 逼近 ， 即 
的 于 2 
J=| LO-7COPw 人 
达到 最 小 。 与 曲线 拟 合 一 样 ， 选 一 组 函数 十 0 二 上 构造 /CD ， 即 令 
JJ=a7iO0+ 十 oO 


代入 (11) 式 ， 求 和 4 使 /达到 极 小 。 利 用 极 值 必要 条 件 可 得 
(1 7 (7) 






























































a) | = 


(7 厂 ) 人 (75 过 03 ( 户 ) ， (10.9) 

















D 
, 硬 三 二 C 
(8 7) 中 解 。 


二 人 
最 简单 的 当然 是 用 多 项 式 吉 近 函数 ， 即 选 i00 =1， 态 CD) = 工 ， 态 (00 三 X ，…。 并 且 如 





这 中 


















































已 

Or Oodaxc=0 ，，， ， 
时 人 作 上 5 本 友 ，(@ 寺 妃 ， 方 程 组 (10.9) 的 系数 矩阵 将 是 对 角 阵 ， 计 算 大 大 简化 。 
满足 这 种 性 质 的 多 项 式 称 正 交 多 项 式 。 
惑 让 得 〈Legendre) 多 项 式 是 在 [1 区 间 上 的 正 交 多 项 式 ， 它 的 表达 式 为 




















1 丰 
石 (oO=1， 有 CD=TNTRCC 一 D 12 
可 以 证 明 
0， 1 尖 7 
1 
BOPCDOa=1 2 
| Ca 1= 7 
27 十 ] 
2K 十] 大 
已， (ZX)= 有 二 1 1 大 三 12,…: 











常用 的 正 交 多 项 式 还 有 第 一 类 切 比 雪夫 〈Chebyshev) 多 项 式 
7 (Z)=cos(1arccoSX)， (xze [一 1 = 0,12…) 
和 拉 盖 尔 (Laguerre) 多 项 式 














也 (z)=e” (xz “e” )，(xe[0,+co), 姑 =0,12…) 
ax” 
JCO=cosxe|- 王 | 有 
例 6 求 2 2. 在 如 =Spanbx ,X 上 中 的 最 佳 平方 逼近 多 项 式 。 
解 ”编写 程序 如 下 : 
SymS 


base=[1,，X^27，X^4]; 
yY1L=base.'xbase 
Y2=Cos (X) xbase. 


zlL=int(yl， -pi/2v pi/2) 





xz2=int(y2，-pi/2， pi/2) 
a 三 TANz2 





Xishul=aqouble (al) 
qidgits (8)，Xxishu2=vpa(al) 


求 得 xishul=0. 9996 -0. 4964 0. 0372， 即 所 求 的 最 佳 平方 通 近 多 项 式 为 
y=0.9996 一 0.4964x 十 0.0372x 


8S 回归 分 析 
回归 分 析 是 处 理 很 难 用 一 种 精确 方法 表示 出 来 的 变量 之 间 关 系 的 一 种 数学 方法 ， 它 是 最 党 
j 的 数理 统计 方法 ， 能 解决 预测 、 控 制 、 生 产 工艺 优化 等 问题 。 它 在 工农 业 生产 和 科学 研究 各 个 
领域 中 均 有 广泛 的 应 用 。 
回归 分 析 一 般 分 为 线性 回归 分 析 和 非 线性 回归 分 析 。 本 节 着 重 介绍 线性 回归 分 析 的 基本 结 
论 及 其 在 Matlab 中 的 相应 命令 。 线 性 回归 分 析 是 两 类 回归 分 析 中 较 简 单 的 一 类 ， 也 是 应 用 较 多 
的 一 类 。 
一 “一 元 线性 回归 分 析 
针对 一 组 〈 二 维 ) 数据 C5， 芒 ) ， 守 = 2 于 〈 其 中 已 互 不 相同 ) ， 其 最 简单 的 数据 拟 
合 形式 为 寻求 直线 了 二 永 二 局 x ， 使 记 十 启 x 在 最 小 二 乘 准 则 下 与 所 有 数据 点 最 为 接近 。 但 
于 随机 观测 误差 的 存在 ， 满 足 上 述 数据 点 的 直线 应 该 是 
7 二 所 TAX+e， (10.10) 
其 中 忆 了 是 准确 的 ，/，, 记 是 两 个 未 知 参数 ，2 是 均值 为 零 的 随机 观测 误差 ， 具 有 不 可 观测 性 ， 
可 以 合理 地 假设 这 种 观测 误差 服从 正 态 分 布 。 于 是 我 们 得 到 一 元 线性 回归 模型 为 
人 
FEe=0,Dse=a?”， 
















































































































































































































































































(10.11) 














5 中 加 未 知 ， 固 定 的 未 知 参数 所 、 记 称 为 回归 系数 ， 自 变量 革 称 为 回归 变量 。 
(10.10) 式 两 边 同 时 取 期 望 得 : 三 = 朱士 记 < ， 称 为 了 对 并 的 回归 直线 方程 。 
在 该 模型 下 ， 第 ?个 观测 值 可 以 看 作 样本 五 三 启 十 万 二 6 《这 些 样本 相互 独立 但 不 同 

分 布 ，z = 了 2 ) 的 实际 抽样 值 ， 即 样本 值 。 

一 元 线性 回归 分 析 的 主要 任务 是 ， 用 实验 值 〈 样 本 值 ) 对 永 、 记 和 G 作 点 估计 ， 对 回归 
系数 万 、 忆 作假 设 检验 ;在 关 = xo 处 对 了 作 预 测 ， 并 对 了 作 区 间 估 计 。 

1 回归 参数 友 、 记 和 0 估计 

假设 有 于 组 独立 观测 值 : (2 六) = 上 2 于 ， 则 由 (10.10) 有 

| 二 二 局 2 十 CE， 
尼 2E; =0,D8 =G 


将 


























































































































其 中 吾 ， = 2 相互 独立 。 记 
=2(0p 0) 了 二 》 (Oy 一 二 记 击 六 
HH 4 




















称 2(po,B) 为 偏离 真实 直线 的 偏差 平方 和 。 由 最 小 一 滋 法 得 到 的 估计 忆 G= 0.D) 称 为 忆 的 最 
小 二 乘 估 计 ， 其 中 


六 人 六 一 1 1 
沁 一 X 一 三 二 志 三 业 2 ~、 工 
记 = 了 一 记 元 必 和 2 了 一 2 》 人 2 
经验 ) 回归 方程 为 


















































了 Y=P+Hxz=7+ACx- 旭 ， ( 作 下 六 
这 样 我 们 得 到 应 公 = 0 册 的 无 偏 估计 应 G= 040) ， 其 中 应 公 = 0 服从 正 态 分 布 


有 N(CD GUI/PT 天 7 y Cr 一 下 有 Na” y Cr 一 避 ? 


且 对 每 组 (5， 芒 ) ， 可 以 求 出 拟 合 值 五 以 及 残 差 疙 一 于 ， 易 知 
》0 -2)=0 
1 。 





























这 说 明 残 差 之 和 为 零 。 
若 记 
2. -2 加 = 乞 志 -应 - 房 zx = 的 人 
则 称 乡 为 残 差 平方 和 或 剩余 平方 和 。 在 模型 (IO.1D 下 ， 易 知 ，ZQ. = O 一 2)G ， 常 记 
0 =0./-2 
称 G。 为 剩余 方差 ( 残 差 的 方差 )， 它 是 G? 的 无 偏 估 计 ， 且 有 如 下 分 布 
] 1 一 2 、， 
0 
9e 分 别 与 思 、 记 独立 ，G。 称 为 剩余 标准 差 ， 显 然 ，9。 越 接近 于 0， 说 明 线性 回归 方程 
(10.12) 越 显著 。 
2 模型 的 假设 、 预 测 、 控 制 
1 回归 方程 的 显著 性 检验 
在 实际 问题 中 ， 因 变量 ”与 自 变 量 Y 之 间 是 否 有 线性 关系 (10.10) 只 是 一 种 假设 ， 在 求 出 回 
归 方程 之 后 ， 还 必须 对 这 种 回归 方程 同 实际 观测 数据 拟 合 的 效果 进行 检验 。 
由 (10.10) 可 知 ，| 忆 越 大 ， 了 随 关 变化 的 趋势 就 越 明 显 ; 反 之，| 尿 | 越 小 ， 了 随 并 变化 的 
趋势 就 越 不 明显 。 特别 当 应 -0 时 ， 则 认为 少 与 基 之 间 不 存在 线性 关系 ， 当 局 关 0 时 ， 则 认为 
与 尉 之 间 有 线性 关系 (10.10)。 因 此 ， 问 题 归 结 为 对 假设 
Zoo:D=0， 已 :Diz=0 
进行 检验 。 假设 如 o : 房 = 0 被 拒绝 ， 则 回归 显著 ， 认 为 了 与 革 之 间 存在 线性 关系 ， 所 求 的 线性 
回归 方程 有 意义 ， 香 则 回归 不 显著 ， 了 与 x 的 关系 不 能 用 一 元 线性 回归 模型 来 描述 ， 所 得 的 回 
归 方 程 也 无 意义 。 此 时 ， 可 能 有 如 下 几 种 情况 : 
GD 并 对 了 没有 显著 影响 ， 此 时 应 丢掉 变量 工 ; 
(ii) 屏 对 儿 有 显著 影响 ， 但 这 种 影响 不 能 用 线性 关系 来 表示 ， 应 该 用 非 线性 回归 
( 诈 ) 除 半 之 外 ， 还 有 其 他 不 可 忽略 的 变量 对 儿 有 显著 影响 ， 从 而 削弱 了 工 对 了 的 影响 。 
此 时 应 用 多 元 线性 回归 模型 。 因 此 ， 在 接受 好 o 的 同时 ， 需 要 进一步 查 明 原因 以 便 分 别处 理 。 




















































































































[本 1 








































































































































































































下 面 介绍 两 种 检验 方法 ， 分 别 是 
(aJF 检 验 法 














对 样本 方差 局 进行 分 解 ， 有 


i=1 






























































上 式 中 的 台 * 是 由 实际 观测 值 没有 落 在 回归 直线 上 引起 的 〈 否 则 为 零 ) ， 习 是 由 回归 直线 引起 
的 。 因 此 ， 忌 越 大 ， 纪 .就 越 小 ， 表 示 了 与 工 的 线性 关系 就 越 显著 ， 否 则 ， 辟 越 小 ， 马 .就 越 大 ， 
表示 】 与 工 的 线性 关系 就 越 不 显著 。 这 样 我 们 就 找到 了 -一 种 判别 回归 直线 拟 合 程度 好 坏 的 方法 ; 
如 果 芝 /1S 接近 于 1， 即 上 /&。 较 大 时 ， 则 对 拟 合 效果 感到 满意 。 

易 知 ， 当 妇 成 立时 ， 有 芝 1G?_Y 四 上 且 纺 与 避 独 立 。 再 由 F 分 布 有 

区 L7 时 (2 一 2)7 过 7 
人 ”1 VO+O， 

其 中 上 称 为 相关 系数 。 对 给 定 的 显著 水 平 & ， 有 置信 水 平 为 1-C 的 临界 值 -c 人 一 所 ， 从 
而 F 检 验 法 的 检验 准则 为 ， 当 瑟 > 和 -db2 一 2 时 ， 拒 绝 刀 。， 否 则 就 接受 如。 。 

(b) t 检 验 法 

当 巡 成 立时 ， 由 T 分 布 的 定义 有 


1 y Cr 52 
太一 i=1 


6 to-2) 







































































































































































防 








此 ， 对 于 给 定 的 显著 水 平 % ， 用 1 统计 量 检验 如。 ， 有 置信 水 平 为 1 一 的 临界 值 
17 一 2 ,从 而 t 检 验 法 的 检验 准则 为: 
当 | 卫 亿 三 wo 人 2 一 2 时 ， 拒 绝 好 。 ， 否 则 就 接受 Lo 。 
2 回归 系数 的 区 间 估计 
广 = 六 Co 一 富 : 
令 和气 ， 则 由 
( 房 -BUGVUn+ 忆 /四 ro-2D，( 房 -BAGVP) xz- 
得 到 在 置信 水 平 1-&% 下 ， 尺 和 记 的 置信 区 间 分 别 为 
[ 房 -aeaCe-26.VLR+ 姑 /， 忆 HH DC-2)G.VL+52 7/ 太 
[应 -opC=-26.7VF 记 +wa-26./WFj， 




































































]， 





和 六 
再 由 丝 /9 _Y (一 2) ， 还 可 得 在 置信 水 平 1 一 w% 下 ，G? 的 置信 区 间 为 
O OO 
[Wireiz02 一 2 ， Ma 一 2]， 




















3) ”预测 与 控制 
当 检验 结果 拒绝 了 习 。: 记 
测 和 控制 。 预 测 就 是 对 固定 的 关 值 预测 相应 的 了 人 
控制 在 制定 的 范围 内 。 
国 
0 并 且 通 过 了 假设 检验 ， 即 检验 结果 拒绝 了 如 o : 尿 =0 。 令 加 
个 国定 值 ， 且 了 > = 办 TDz+Eo， E0_N(0a”) ,假设 yo 区 相互 独立 ， 

















= 0 ， 接 下 来 的 问题 是 如 何 利用 回归 方程 了 = io + 忆 x 进 行 预 
控制 就 是 通过 控制 的 值 ， 以 便 把 了 的 值 
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表示 并 的 某 

na 

加 的 预测 值 为 加 的 回归 值 下 三 启 十 启 zo 。 它 是 瑟 o 的 无 偏 估计 ， 即 
E = ELB+ 记 = 证 


纷 定 显著 水 平 & ，Jm 的 轻信 水 平 为 1- C 的 预测 区 间 为 [0 一 9(xo) ， 和 0 十 (xzo)]， 其 中 
ou)=GNO-2)VJI+LAAT+(O 一 宫 “/ 











































































































上 式 可 知 ， 剩 余 标准 差 C 越 小 ， 预 测 区 间 越 小 ， 预 测 值 越 精确 ， 对 于 给 定 的 样本 观测 值 和 
置信 水 平 而 言 ，”。 越 靠近 三 时 ， 预 测 精 度 就 越 高 。 


由 于 2 的 任意 性 ， 夹 在 两 曲线 
六 (xo ) = 六 -Or) ， 六 (zxo) = 六 +0O(xo) 



















































































之 间 的 部 分 就 是 了 = 上 + 万 x+e 的 置信 水 平 为 1- C 的 预测 带 。 特 别 ， 当 于 很 大 且 工 在 天 附近 
取 值 时 ， 有 


























证人 2 二 古 VJL+LP+(Om 一 祝 27 太 = 下 
其 中 乓 -ea 是 置信 水 平 为 1 CX 的 标准 正 态 分 布 的 临界 值 。 此 时 ， (CD = ce-a'a ， 于 是 了 的 
置信 水 平 为 1 一 C 的 预测 区 间 近 似 为 [” 一 Core ， ”十 Core2]。 在 这 种 情况 下 ， 预 测 带 是 
平行 于 回归 直线 的 两 直线 之 间 的 部 分 ， 预 测 工作 得 到 了 很 大 的 简化 。 例 如 ， 在 应 用 中 可 简单 地 
取 CQ =0.03 ， 则 允 -ev2 一 %oo7s 一 1.90 ;， 取 CQ = 0.01， 则 史 -ev2 二 %ooos 一 2.38 .从 而 得 到 署 信 
水 平 为 95% 与 99% 的 预测 区 间 分 别 近似 为 [了 一 26C。， 了 二 1 3060] 和 [了 “580。， 
了 +2.580.。]。 

人) 控制 

车 要 了 “二房 十 厅 x+E 的 值 以 1-- C 的 概率 落 在 指定 区 间 (c,d) 之 内 ， 变 量 x 应 控制 在 什 
入 范围 内 的 问题 就 是 所 谓 的 控制 问题 。 它 是 预测 问题 的 反问 题 。 

由 预测 的 讨论 知 ， 夹 在 两 曲线 五 (xo) 二 加 一 9(xo)，2(Cxo) 三 二 9(Cxo) 之 间 的 部 分 就 
是 了 三 启 十 记 x+8 的 置信 水 平 为 1- C 的 预测 带 ， 故 若 要 了 的 观测 值 以 1- C 的 概率 落 在 区 间 
0 

-OOD)>c， 了 +OoCD)<d 
这 要 求 @& 一 c>26(xz) ， 若 方程 了 -6(COo = C ， 了 +6(Oo = Q 分 别 有 解 纪 迟 ， 则 (Q,D) 就 是 所 求 





































































































































































































的 工 的 控制 区 间 。 

一 般 来 说 ， 要 解 出 @ 和 睛 很 复杂 。 若 样本 容量 很 大 ， 即 丑 很 大 且 习 在 亏 附 近 取 值 时 ， 一 个 
简化 的 方法 就 是 利用 近似 预测 区 间 [了 一 人 eu-eia ， 了 十 Gewi-ei?] 来 进行 控制 。 

注 ， 一 元 线性 回归 是 后 面 介 绍 的 多 于 线性 回归 的 特例 ， 其 Matlab 实 现 见 本 节 第 五 部 分 。 

一 “可 线性 化 的 一 元 非 线性 回归 《曲线 回归 ) 

在 工程 技术 中 ， 自 变量 x 与 因 变 量 / 之 间 有 时 呈现 出 非 线性 〈 或 曲线 ) 关系 ， 这 是 通常 出 
现 两 种 情况 :一 种 是 呈现 多 项 式 的 关系 ， 这 种 情况 通过 变量 蔡 换 可 化 为 多 元 线性 回归 问题 给 巴 
解决 ， 另 一 种 是 呈现 出 其 它 非 线性 关系 ， 通 过 变量 替换 可 化 为 一 元 线性 回归 问题 给 予 解决 。 
已 知 变量 x 与 少 的 一 组 〈 二 维 ) 观测 数据 @5 ;六 ) ，= 0 于 《其 中 交互 不 相同 ) ， 
通常 做 出 二 维 平面 上 的 散 点 图 ， 着 散 点 的 分 布 不 接近 于 一 条 直线 ， 通 常 根据 散 点 图 的 类 型 匹配 
相应 的 拟 合 曲线 。 若 匹配 曲线 为 多 项 式 

JI=D+DXZT TD (7 <7) ， 


这 是 线性 回归 问题 ， 即 关于 回归 系数 永 (= 0 …, 到 ) 是 线性 的 。 此 时 只 需 令 


人 5 
放大， 



















































































































































































































































































































































































尼 
杞 








归 问 题 化 为 多 元 线性 回归 问题 ， 其 解决 方法 详 见 下 面 介绍 的 多 元 线性 回归 问题 〈 另 外 ， 多 
项 式 回 归 问 题 的 实现 在 Matlab 中 有 相应 的 命令 ， 有 具体 见 多 元 线性 回归 部 分 》 。 关 匹 配 曲线 〈 经 


验 公式 ) 为 含 参量 4 的 非 线性 曲线 ， 采 用 的 办 法 是 通过 变量 蔡 换 把 非 线性 回归 化 为 线性 回归 。 









































































































































通常 匹配 的 含 参量 人 之 的 非 线 性 曲线 有 以 下 六 类 ， 具 体 的 蔡 换 方法 如 下 
人 
1 双 曲 线 具 
1 1 
1 三 一 三 一 2 三 一 
作 变 量 蔡 换 。 X， 了， 得 Y= 4+bu 。 由 观测 值 CD) = 0 2 按 厂 ， 
] 
RE 

















六 算出 人 六) ， = 0 。 对 妊 与 ?利用 前 面 的 经 验 回归 直线 公式 ， 计 算出 参数 估计 
值 记 、 启 ， 则 有 回归 方程 
] 局 

X 


















































2 过 函数 曲线 了 = 0 ， 其 中 Z> 04>0. 
两 边 取 常 用 对 数 ，]8) 二 可 4 + 0 井 苹 ， 再 作 代 换 世 一 加 X，Y 二 区》， 作 二 塌 4 ， 则 寡 函 
数 田 线 方程 就 变 成 直线 方程 





















































7V=4+pDu ， 
由 观测 值 C2， 疙 ) ， 1 = 0 于 按 帮 三 吉 疙 ， 疙 二] 六 算出 GoV) ， 0 对 
4 与 利用 本 节 前 面 的 经 验 回归 直线 公式 ， 计 算出 参数 估计 值 4、 ， 又 由 6@=e“ 得 到 4 。 因 此 
































有 回归 方程 

















7 了 =bxz” 
3 指数 函数 曲线 一 ae" ， 其 中 > 0 。 
两 边 取 常用 对 数 ， 井 三 4+Ax ， 再 作 代 换 & = 并 ,7 二 1] ，4 二 井 4 ， 则 指数 函数 
曲线 方程 就 变 成 直线 方程 























7V=4+pDu ， 
观测 值 Co] ) ， 1 0 天 按 臣 二 半 ， 疙 三 井 世 算出 (7) ， 人 对 1 
与 利用 本 节 前 面 的 经 验 回 归 直 线 公式 ， 计 算出 参数 估计 值 4 2 ， 又 由 @= e“ 得 到 .因此 有 
回归 方程 





































































































y=be” 
4 倒 指 数 函数 曲线 了 = ae “， 其 中 &> 0 。 
D 1 
1]1g7y =1lga+ 一 4U 三 一 
两 边 取 常 用 对 数 ; x ， 再 作 代 换 。 x，Y=187，4= 雹 4 ， 则 倒 指 数 函 
数 曲线 方程 就 变 成 直线 方程 
7V=4+pDx ， 




















由 观测 值 Cp 攻 ) = 0 和 光 按 格 一 ] 攻 ， 玫 二 本 力 算 出 他 01… 光 。 对 
2 放 4 


入 


4 与 Y 利 用 本 节 前 面 的 经 验 回归 直线 公式 ， 计 算出 参数 估计 值 人 .2 ， 又 由 4=e 得 到 4 .因此 
有 回归 方程 



























































六 五 7 天 
少 二 CE 


5 对 数 函 数 曲 线 了 = 4+lgxY ， 其 中 xz> 0。 
作 代 换 & = 区 x ，Y 三 隐 ， 则 对 数 函数 曲线 方程 就 变 成 直线 方程 


Vv=Q+pu 


























观测 值 00， 蕊 ) ， 1 0 天 按 大 三 雹 六 ， 呈 一 沁 算出 亿 o7) ， 1 = 0 。 对 
与 利用 本 节 前 面 的 经 验 回归 直线 公式 ， 计 算出 参数 估计 值 4 、& 。 由 此 得 到 回归 方程 


















































































































































=0+Dlgx 
1 
儿 二 三 
6 S 型 函数 曲线  Q+pe 。 
1 攻 1 
一 =4+pDe 本 
上 述 方程 等 价 于 了 ， 作 代 换 &= ee 一， 卫 ， 则 S 型 函数 曲线 方程 就 变 成 直线 
方程 
7v=Q+pDU ， 
由 观测 值 C? 了 = 0 天 算出 Cs) 04 于。 对 刀 与 利用 本 节 前 面 的 经 
验 回归 直线 公式 ， 计 算出 参数 估计 值 、 。 因 此 有 回归 方程 
1 
+ben 





注 : 对 于 非 线性 回归 问题 的 Matlab 实 现 问题 ， 一 种 方法 是 化 为 相应 的 线性 模型 实现 ， 另 一 



































































































































































































































种 方法 是 直接 应 用 Matlab 中 相应 的 命令 ， 其 结果 是 一 致 的 。 详 见 本 节 第 五 部 分 。 
三 “多 元 线性 回归 分 析 
一 般 地 ， 在 实际 问题 中 影响 应 变量 ] 的 自 变量 往往 不 止 一 个 ， 不 妨 设 有 大 个 为 ，722 ， 
…，x 。 通 过 观测 得 到 一 组 (大 十 1 维 ) 相互 独立 的 试验 观测 数据 Ca 六 07 了 7 ， 
7 2 …， 其 中 站 > 大 +]。 假 设 变量 与 变量 六 ， 浆 ，…， 交 之 间 有 线性 关系 ; 
=D+DDT… 十 x+e， (10.13) 
其 ， 人 0，De =a-” > 0。 则 观测 数据 满足 
= [+ 万 2 十 … 二 二 E， ) 也 三 12,…: (10.14) 
人 
JJ 1 2 1 D El 
= 区 克 二 的 12 68= 
JJ， ] ln 光 D CE， 
则 (10.14) 可 简写 为 
了 =X15+ae， (10.15) 
其 中 无为 已 知 的 有 X (十 芒 矩阵 ， 称 为 回归 设计 矩阵 或 资料 矩阵 ， 了 是 于 维 观察 值 列 向 量 ， 
太 为 上 + 工 维 未 知 的 列 向 量 ，2 是 满足 
瑟 E=0)， 
COF(e,e)=a27 
的 球 维 随机 列 向 量 ， 其 中 0 为 二 维 零 向 量 ，G? 是 未 知 参数 ， 2 ，7=12 
为 严 阶 单位 矩阵 ， 即 对 随机 误差 2 ，…，& 作 无 偏 、 等 方差 与 互 不 相关 的 假定 。 一 般 称 
了 =XO+2， 
尼 E =0)， 
1 
COF es=0 1 0010) 
为 大 线性 回归 模型 《高 斯 一 马尔 科 天 线性 模型 ) ， 并 简 记 为 (to 
对 (10.16) 取 数学 期 望 得 到 


称 为 线性 回归 方程 。 


三 记 十 启 辣 十 二 Dr ， 











对 线性 模型 (ZXB,a 闷 ) ， 所 要 考虑 的 主要 问题 是 ; 
(Gi 用 实验 观测 数据 对 未 知 参数 记 ，…， 
































量 和 自 变量 阅 ，…， 之 间 的 线性 关系 ; 
(ii) 在 冯 一 ao 丸 二 <60 处 对 上 的 值 作 预 测 和 控 
是 假设 导 > 上 +1 。 





友和 Ga? 做 点 估计 和 假设 检验 ， 从 而 建立 因 


制 ， 





变 











并 对 了 作 区 间 估 


计 。 本 部 分 总 











1 未 知 参数 房 和 Ga? 估计 
先 用 最 小 二 乘法 求 启 ，.…， 计 的 估计 量 。 做 误差 平方 和 


Os》 OU -0 一 Dr 0 


























首 





























， 永 作为 房 ，…， 永 的 点 估计 。 由 微分 学 





寻求 使 得 达 到 最 小 值 的 ， 值 的 方法 


知 ， 上 只 需求 解法 方程 组 《或 称 正 规 方程 组 ) 










































































































































































































































































和 1=X 了 ， 
其 中 瑟 、 了 同 前 ， 太 = (5 房 ，… 记 ) 。 解 的 估计 值 太 = (人 台 ) 王子 ， 这 样 得 到 的 永 带 
= 记 +DD+… DB 
in 同时 有 结论 
) 启 服 从 正太 分 布 ， 尼 服从 丰 十 1 维 正 态 分 布 
Ci) 及 是 记 的 无 偏 估计 且 是 最 优 线 性 无 伪 估 计量 ， 协 方差 矩 阵 为 C_(X 只 ) 。 
0 2=G-Gh)TGT-GD) 共和 G=[haez] 出 有 EOJ=O-E-Da2 
且 
池 . 
ao2 _X (2 大-]) 
O O 一 他 ? 
从 而 1-E-L = ， 即 1-E-1 为 的 无 偏 估计 量 。G 是 剩余 方差 ，G 是 剩余 标 
准 差 
三 (六 一 羽 ” 
Civ) 对 了 的 样本 方差 气 ， 有 
人 
So4U 地 忆 
其 中 包 是 残 差 平方 和 ， 反 映 随 机 误差 对 了 的 影响 ，C 是 由 回归 方程 引起 的 ， 称 为 回归 平方 和 ， 
反映 自 变 量 对 了 的 影响 。 
2 多 元 线性 回归 中 的 假设 检验 
在 实际 问题 中 ， 往 往事 先 不 知道 或 不 能 确定 随机 变量 》 与 自 变量 国 ，…， ?之 间 确 有 线 
性 关系 。 因 而 (10.13) 往 往 是 一 种 假设 ， 因 此 在 求 出 线性 回归 方程 之 后 ， 还 必须 对 求 出 的 线性 回 
月 广 各 同 实际 观测 数据 所 全 类 似 于 一 元 线性 回归 ， 可 提出 以 下 原 假设 
已 0 : 0 = =…=A=0 
L7 
当 人 既成 立时 ， 着 令 。 【C++ 丝 ( 称 尺 为 相关 系数 ) ， 则 有 
Un-1 
Op-kt-D 人 1-R2 Fa-E-D 








当 显著 水 平 忆 给 定 后 ， 由 民 一 人 > 人 -人 7 了 | 人 成语 
1-w= 串 严 < (La 天 -TD| 丽 成 立 | 


， 即 














知 F 检 验 法 的 检验 规则 为 ; 

如 果 轿 > 三 -et 一 大 一 芒 ， 则 拒绝 好 。， 认 为 因 变量 与 自 变 量 刁 ，…， 玖 之 间 的 线 
性 关系 显著 ;否则 ， 认 为 了 与 嫉 ，…， 关 之 间 的 线性 关系 不 显著 。 

需要 注意 的 是 ，) 了 与 ，…， 之 间 的 线性 关系 不 显著 ， 可 能 出 现 几 种 情况 ， 如 了 于 其 
中 某 些 自 变 量 无 关系 ， 可 以 去 掉 这 些 自 变量 ;了 与 为 ，…， 人 反之 间 的 存在 非 线 性 关系 ;还 有 
其 它 变量 与 了 有 关系 等 。 当 然 还 有 其 它 检 验方 法 。 

3 多 元 线性 回归 中 的 预测 

Gb 点 预测 

当 我 们 求 出 回归 方程 


























































































































》》=D+DT TD 
并 经 过 检验 后 ， 对 于 给 定 的 自 变量 的 值 石 一 Xio ，…， 一 20， 自然 用 加 一 记 十 语 mno 
十 十 应 xao 来 预测 Jo=D+A2o+… Dixio+e 称 加 为 加 的 点 预测 。 实际 上 ， J 加 是 
己 。 的 无 偏 估计 量 〈 因 为 《Oo mo) 一 0 ) 
(Gi) 区 间 预 测 
由 于 及 Na(Da (CT TD) 所 以 】- N(GP,Ga ( 莹 GD) 当 G 已 知 时 ， 由 
于 纪 与 记 相 五 独立 ， 从 而 2 与 了 = 人 


) 与 马 相 互 独立 。 从 而 少 一 少 与 马 相 互 独立 ， 且 有 
JJ 一 了 Na +CGCX TD GD) 
















































































令 | 了 一 少 医 开 ， 即 了 的 置信 区 间 为 少 一 天 , 少 + 天 ) ， 天 由 置信 度 1 一 C 确定 ， 则 当 : 


司 
踢 



























































台 
水 平 % 给 定 后 ， 根据 G ~X (下 
六 二 2 
aWJI+GCrT7Z)TG7) 0 t-D 











从 而 
关 二 靖 全 三 大 1 风 人 学 和: 














此 当 显 著 水 平 C% 给 定 后 ， 了 的 置信 区 间 为 消 一 到 , 销 十 天) ， 其 中 天 由 上 式 确定 。 
注 ， 多 元 线性 回归 的 Matlab 实 现 详 见 本 节 第 五 部 分 。 

四 “逐步 线性 回归 分 析 
从 多 元 线性 回归 分 析 中 我 们 知道 ， 采 用 的 自 变量 越 多 ， 则 回归 平方 和 越 大 ， 残 差 平方 和 越 
小 。 然 而 ， 采 用 较 多 的 变量 来 拟 合 回归 方程 ， 得 到 的 方程 稳定 性 差 ， 每 个 自 变量 的 区 间 误差 的 
识 累 将 影响 总 体 误差 ， 用 这 样 建立 起 来 的 回归 方程 作 预 测 的 可 靠 性 差 、 精 度 低 。 另 一 方面 ， 如 果 
采用 了 对 因 变 量 影响 小 的 自 变量 而 遗漏 了 重要 变量 ， 可 导致 佑 计量 产生 偏 倚 和 不 一 致 性 。 因 而 
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岗 











希望 得 到 最 优 的 回归 方程 。 
逐步 线性 回归 分 析 方 法 就 是 一 种 自动 从 大 量 可 供 选 择 的 变量 中 选择 那些 对 建立 回归 方程 比 
较 重 要 的 变量 的 方法 ， 它 是 在 多 元 线性 回归 基础 上 派生 的 一 种 算法 技巧 ， 详 可 参阅 相应 的 文献 。 
其 基本 思路 为 : 从 一 个 自 变量 开始 ， 视 自 变量 对 ] 作用 的 显著 程度 ， 从 大 到 小 依次 逐个 引入 回 
归 方 程 。 当 引入 的 自 变 量 由 于 后 面 自 变 量 的 引入 而 变 得 不 显著 时 ， 要 将 其 剔除 掉 。 引 入 一 个 自 变 
量 或 从 回归 方程 中 剔除 一 个 自 变量 ， 为 逐步 回归 的 一 步 。 对 于 每 一 步 ， 都 要 进行 》 值 检验 ， 以 
保 每 次 引入 新 的 显著 性 变量 前 回归 方程 中 只 包含 对 作用 显著 的 变量 。 这 个 过 程 反复 进行 ， 
直至 即 无 不 显著 的 变量 从 回归 方程 中 史 除 ， 又 无 显著 变量 可 引入 回归 方程 止 。 
使 用 逐步 线性 回归 时 要 注意 : 要 适当 选择 引入 变量 的 显著 性 水 平和 剔除 变量 的 显著 性 水 平 ; 
应 尽量 选择 那些 相互 独立 性 强 的 变量 。 
五 回归 分 析 的 Matlab 实 现 
Matlab 统 计 工 具 箱 中 提供 了 一 些 回 归 分 析 的 命令 ， 现 介绍 如 下 。 
1 多 元 线性 回归 
多 元 线性 回归 的 命令 是 regress， 此 命令 也 可 用 于 一 元 线性 回归 。 其 格式 为 : 
(1) 确定 回归 系数 的 点 估计 ， 用 命令 : b=regress(Y，X)。 
《2) 求 回 归 系 数 的 点 估计 和 区 间 人 估计， 并 检验 回归 模型 ， 用 命令 : 
[b，bint，Tr，Trint，stats]=regress(Y，X，alpha)。 
3) 画 出 残 差 及 其 置信 区 间 ， 用 命令 : 
Tcoplot (CT，Trint) 。 
在 上 述 命令 中 ， 各 符号 的 含义 为 : 
GD 2 三 及，Y，X 的 定义 同 本 部 分 前 面 所 述 。 对 一 元 线性 回归 ， 在 加 三 记 ，Y，X 中 取 大 
=]1 即 可 ; 
(ii) alpha 为 显著 性 水 平 〈 缺 省 时 为 0.05) ; 
《iii) bint 为 回归 系数 的 区 间 估 计 ; 
Gv) r 与 rint 分 别 为 残 差 及 其 置信 区 间 ; 
Cv) stats 是 用 于 检验 回归 模型 的 统计 量 ， 有 三 个 数值 ， 第 一 个 是 及 ， 第 二 个 是 F 值 ， 第 
三 个 是 与 FE 对 应 的 概率 P。 其 中 尽 ” 与 F 定 义 同 前 ， 值 越 大 ， 说 明 回 归 方程 越 显著 ， 己 < % 时 拒绝 
生 。， 回 归 模型 成 立 。 
例 7 合金 的 强度 /与 其 中 的 碳 含量 有 比较 密切 的 关系 ， 今 从 生产 中 收集 了 一 批 数据 如 下 
攻 。 试 先 拟 合 一 个 函数 Jo0 ， 归 分 析 对 它 进 行 检验 。 
X 0.10 0.11 0.12 0.13 0.14 0.13 0.16 0.17 0.18 
y 42.0 41.3 45.0 43.3 43.0 47.3 49.0 33.0 50.0 
解 ” 先 画 出 散 点 图 : 
X=0.10:0.01:0.18; 
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y=[42.0,41.3,45.0,45.3,43.0,47.3,49.0,35.0,50.0]; 


polt(xcy, 十 ) 


可 知 y 与 x 大 致 为 线性 关系 。 
设 回归 模型 为 ?= 十 局， jregress 和 rcoplot 编 程 如 下 : 





clc,clear 


xl1=[0.10:0.01:0.18] ; 
y=-[42.0,41.5,45.0,45.5,45.0,47.5,49.0,55.0.50.0] ; 


X=[ones(9,1),X1]; 



































[b,bintrrint,stats]=regress(YX); 


b,bintstats,rcopolt(Print) 


得 到 
b=27.4722 


bint=18.6851 
73.7735 


stats=0.7985 
即 加 -274722， 太 -137.5000， 议 的 置信 区 间 是 [18.6851， 36.2594]， 忆 的 置信 区 间 是 


137.3000 


306.2594 
199.2245 


27.7469 ”0.0012 


























[75.7755，199.2245]; 及 "=-0.7985， 下 =27.7469， 刀 三 0.0012。 
可 知 所 设 回 归 模 型 成 立 。 
观察 命令 rcoplot(prinb 所 画 的 残 差 分 布 ， 除 第 8 个 数据 外 其 余 残 差 的 置信 区 间 均 包含 零点 ， 





第 8 个 点 应 视 为 蜡 常 点 ， 将 其 别 除 后 


b=30.7280 


























新 计算 ， 可 得 


由 











109.39835 


bint=26.2805 


76.9014 
stats=0.9188 





























35.2834 
141.8955 
67.8534 0.0002 
应 该 用 修改 后 的 这 个 结果 。 
































例 8 某 厂 生产 的 一 种 电器 的 销售 量 ” 与 竞争 对 手 的 价格 x1 和 本 厂 的 价格 x2 有 关 。 下 表 是 该 














厂商 品 在 10 个 城 了 














的 销售 记录 。 试 根据 这 些 数 据 建 立 y 与 x1 和 x2 的 关系 式 ， 对 得 到 的 模型 和 系数 






































进行 检验 。 若 某 市 本 三 产品 售 价 160《〈 元 ) ， 涡 争 对 手 售 价 170《〈 元 ) ， 预 测 商品 在 该 市 的 销售 


上 了 


虽 
曲 
o 





xl 元 ”120 140 190 130 155 175 125 145 180 150 


X2 元 
Y 个 


解 分 别 画 


100 110 90 130 210 150 250 270 300 2530 


102 100 120 77 46 93 26 09 65 85 











之 间 的 关系 则 


任 以 胡 























Hy 关 于 x1 和 y 关 于 x2 的 散 点 图 ， 可 以 看 出 y 与 X2 有 较 明 显 的 线性 关系 ， 而 y 与 xl 

















定 ， 我 们 将 作 几 种 演 试 ， 用 统计 分 析 决 定 优 劣 。 











设 回 归 模 型 为 
7y= 有 +Rz+Dm 
编写 如 下 程序 ; 
Xl=[120 140 190 130 153S3 175 1253 143 180 1530] | 
Xx2=[100 110 90 130 210 130 230 270 300 250] ; 
=-[102 100 120 77 46 93 26 69 65 85] ; 
X=[ones(10,1),x1,x2]; 
[b,bintprint,stats]=regresS(Y,X); 
b,bint,stats 
得 到 
b=00.3176 0.4139 -0.2698 
bint=-32.30600 163.5$411 
-0.2018 ”1.0290 
-0.4611 -0.0785 
Stats=0.0327 “0.3786 0.0247 
可 以 看 出 结果 不 是 太 好 ，p=0.0247， 取 & =0.05 时 所 设 回 归 模 型 可 用 ， 但 取 C& =0.01 时 所 设 
回归 模型 不 能 用 ， 尺 -0.6527 较 小 ; 所 ， A 的 置信 区 间 包 含 了 零点 。 下 面 将 试图 用 x1，x2 的 
二 次 函数 改进 它 。 
2 多 项 式 回归 
(1) 一 元 多 项 式 的 回归 和 预测 
一 元 多 项 式 的 回归 和 预测 可 用 命令 polyfit 或 polytool 和 polyval 或 polyconf 来 实现 。 殿 
如 下 : 
令 民 二 Go 二 On) 忆 =(002 GoH) 是 多 项 式 
?07 二 Oo ”十 二 onx 十 gob 的 系数 ，8 是 一 个 矩阵 ， 用 来 估计 预测 误差。 
归 可 以 用 命令 [PS]=polyfit(x,ym) 或 polytool(Cxyam) 实 现 ， 其 中 [PS]=polyfit(x,ym) 是 确定 
多 项 式 系数 的 命令 ;命令 polytool(x,yim) 命 令 产生 一 个 交互 式 的 画面 ， 在 画面 中 绿色 曲线 为 拟 
合 曲线 ， 它 两 侧 的 红线 是 了 的 置信 区 间 。 可 以 用 鼠标 移动 图 中 的 十 字 线 来 改变 图 下 方 的 x 值 ， 也 
可 以 在 窗口 内 输入 ， 左 边 就 给 出 y 的 预测 值 与 置信 区 间 。 通 过 左下 方 的 Export 下 拉 式 菜单 ， 可 以 
输出 回归 系数 等 。 
贡 测 和 预测 误差 估计 的 命令 为 polyval 或 者 polyconf， 其 中 Y=polyval(p,x) 求 polyfit 所 得 的 回 
归 多 项 式 在 x 处 的 预测 值 Y;[YDELTA]=polyconftp,x,S,alpha) 求 polyfit 所 得 的 回归 多 项 式 在 x 处 
的 预测 值 Y 及 预测 值 的 显著 性 为 1-alpha 的 置信 区 间 Y 土 DELTA;，alpha 缺 省 时 为 0.05。 
一 元 多 项 式 回 归 也 可 化 为 多 元 线性 回归 来 解 。 
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例 9 观测 物体 降落 的 距离 ”与 时 间 关 的 关系 ， 得 到 数据 如 下 表 ， 求 了 关于 头 的 回归 方程 
了 = QTDxr+cx” 





x 1/30 2/30 3/30 4/30 5/30 6/30 7/30 
y 11.86 15.67 20.60 26.67 33.71 41.93 51.13 
x 8/30 9/30 10/30 11/30 12/30 13/30 14/30 
y 61.49 72.90 85.44 99.08 113.77 129.54 146.48 
解 方法 〈 一 ) “用 一 元 多 项 式 回 归 ， 编 写 程序 如 下 : 
x=1/30:1/30:14/30; 
y=[11.86 15.67 20.60 26.67 33.71 41.93 51.1361.49 72.90 85.44 99.08 
113.77 129.54 146.48]; 
[p,S]=polyfit(x,y,2); 
得 到 p= 
489.2946 65.8896 9.1329 
即 a=9.1329 b=65.8896 c=489.2846。 
方法 〈 二 ) ”化 为 多 元 线性 回归 ， 其 程序 为 : 
x=1/30:1/30:14/30; 
y=[11.86 15.67 20.60 26.67 33.71 41.93 51.1361.49 72.90 85.44 99.08 

















113.77 129.54 146.48]; 
T=[ones(14,1),tt.^2] 
[b,bintprint,stats]=regress(y,T); 
b,stats 
得 到 结果 : b=9.1329 65.8896 489.2946 
Stats 王 1.0e+T007 
0.0000 1.0378 0 
可 以 看 出 ， 两 种 方法 的 出 的 结果 是 一 样 的 。 
(2) 多 元 二 项 式 回 归 
二 项 式 回 归 可 用 命令 : rstool(xymodelalpha)。 其 中 ， 输 入 数据 x\y 分 别 为 姻 X 殉 和 窃 阵 和 
人 alpha 为 显著 性 水 平 〈 缺 省 时 为 0.053) ; model 由 下 列 4 个 模型 中 选择 1 个 “〈 用 字符 串 


和 输入 ， 缺 省 时 为 线性 模型 ) : 
者 

































































yy=D+DT DO + 7 
purequadratic 〈 纯 二 次 ) : 
下 二 DA +… 十 DBx Xm 十 re。 


interaction (交叉 ) : JS 你 


= 克 +JN+ 二 了 人 

quadratic〈 完 全 二 次 ) : 1S7ASm 

命令 rstool 产 生 一 个 交互 式 画 面 ， 画 面 中 有 m 个 图 形 ， 这 m 个 图 形 分 别 给 出 了 一个 独立 变量 
xi《〈 另 m-1 个 变量 取 固 定 值 》 与 y 的 拟 合 曲线 ， 以 及 y 的 置信 区 间 。 可 以 通过 键入 不 同 的 xi 的 值 来 
获得 相应 的 y 值 。 

图 的 左下 方 有 两 个 下 拉 式 菜单 ， 一 个 荣 单 Export 用 以 癌 Matlab 工 作 区 传送 数据 ， 包 括 
beta〈 回 归 系 数 ) 、rmse〈 剩 余 标准 差 ) 、residuals 〈 残 差 ) 。 另 一 个 菜单 model 用 以 在 上 述 4 个 模 
型 中 选择 。 可 以 分 别 选 4 个 模型 ， 并 比较 它们 的 剩余 标准 差 ， 其 中 最 接近 于 0 的 模型 是 最 好 的 。 

我 们 再 作 一 遍 例 8 商品 销 售 量 与 价格 问题 ， 选 择 纯 二 次 模型 ， 即 

=D+Dm+Prt+DOTD ， 
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编程 如 下 : 

xl=[120 140 190 130 155 175 125 145 180 150] ; 
x2=[100 110 90 150 210 150 250 270 300 250] ; 
y-[102 100 120 77 46 93 26 69 65 85] ; 

















X=[xXl x2]; 
Tstool(X,y， puredquadratic ) 
得 到 一 个 交互 式 画 面 ， 给 出 两 幅 图 形 。 左 边 图 形 是 xl 固定 时 的 曲线 yCxJ) 及 其 置信 区 间 ， 右 边 图 






























































形 是 x2 固 定时 的 曲线 y(x2) 及 其 置信 区 间 。 用 鼠标 移动 图 中 的 十 字 线 ， 或 在 图 下 方 窗口 内 输入 ， 
可 改变 xl、x2.。 画 面 左 边 给 出 y 的 预测 值 即 其 置信 区 间 ， 用 这 种 画面 可 以 回答 例 8 提 出 的 “ 若 谋 市 
本 厂 产品 售 价 160〈 元 ) ， 竞 争 对 手 售 价 170〈 元 ) ， 预 测 商品 在 该 市 的 销售 量 ” 问 题 。 

在 画面 左下 方 的 下 拉 式 菜单 Export 中 选择 “all”， 则 beta、rmse 和 residuals 都 传送 到 Matlab 
工作 区 中 。 在 Matlab 工 作 区 中 输入 命令 : 






























































































































































beta,rmse 
得 到 结果 : beta=-312.$871 7.2701 -1.7337 -0.0228 0.0037 
Tmse=10.0430 





如 果 在 另 一 荣 单 model 选 择 其 它 多 元 二 项 式 模型 ， 比 较 它 们 的 剩余 标准 差 就 会 发 现 ， 本 例 
的 所 选 模型 的 rmse=16.6436 最 小 。 

注 : 本 例 中 的 模型 亦 可 化 为 多 元 线性 回归 来 做 。 请 读者 自己 编程 并 比较 结果 。 

3 非 线性 回归 
非 线性 回归 可 用 命令 nlinfitnlintooLnlparcinlpredci 来 实现 。 命 令 格式 如 下 : 
回归 : 回归 可 用 命令 [betarJ]=nlinfit(x,ymodelbeta0) 或 者 nlintool(x,ymodelbeta0,alpha) 
来 实现 。 其 中 命令 [betar 相 =nlinfit(x,ymodelLbeta0) 的 作用 为 确定 回归 系数 ， 而 命令 
nlintool(x,ymodeLbeta0,alpha) 产 生 一 个 交互 式 的 画面 ， 画 面 中 有 拟 合 曲线 和 y 的 置信 区 间 。 通 过 
左下 方 的 Export 下 拉 式 染 单 ， 可 以 输出 回归 系数 等 。 






























































































































































这 里 的 输入 数据 x、y 分 别 为 姑 关 玉 算 阵 和 n 维 列 向 量 ， 
model 是 事先 用 m- 文 从 











T《〈 残 差 ) 





测 


记 >v 



































定义 的 非 线性 函 
\、J (Jacobian 狼 
预测 和 预测 误差 估 


E 阵 ) 是 估 



























































某 些 非 线性 回归 也 可 化 为 多 元 线性 
例 10 在 看 
式 为 
3 0 
其 中 订 ，.…， 永 式 未 知 系数 ，2s72*z 是 
是 反应 速度 。 今 测 的 一 
访 的 参考 值 为 《0.1.0.05.0.02.12) 。 
序号 反应 速度 y 
1 8.55 470 
2 3.79 ”285 
3 4.82 ”470 
4 0.02 ”470 
5$ 2.7$ 470 
6 14.39 ”100 
7 2.54 100 
8 4.353 470 
9 13.00 100 
10 8.50 100 
11 0.05 100 
12 11.32 285 
13 3.12 285 














解 首先， 以 回归 系数 和 自 变量 为 输入 变量 











huaxue.m: 





function yhatF=huaxue(beta.x); 
yphat=(beta(4)#*x(2)-X(3)/beta(3))./(1+beta(])*x(]1)+beta(2)*#x(2)+beta(3)*#x(3)); 





然后 ， 用 nlinfit 计 算 回 归 系 数 ， 月 
即 其 置信 区 间 ， 编 程 如 下 ; 























clc,clear 

X0=[] 8.S5 
2 3.79 
3 4.82 
4 0.02 





470 
285 
470 
470 





300 
80 
300 
80 





一 3 /0 
1+ 启 辣 +Dxa TD 





组 数据 如 下 表 ， 试 | 


300 
80 
300 
80 
80 
190 
80 
190 
300 
300 
80 
300 
190 








10 

10 
120 

120 


数 ; beta0 是 






































对 一 元 非 线 性 回 
回归 系数 的 初 值 .Beta 是 估计 出 的 回归 系数 ， 
计 预 测 误差 需要 的 数据 。alpha 为 显著 性 水 平 ， 缺 省 时 为 0.05。 
计 : 预测 和 预测 误差 估计 的 命令 格式 为 

betaci=nlparci(betapJJ) 其 用 途 为 记 算 回归 系数 的 置信 区 间 ; 

[YDELTA]=nlpredci(modelx,betarJ) 其 
直 Y 及 预测 值 的 显著 性 为 1-alpha 的 置信 区 间 Y 士 DELTA; 
回归 来 解 。 














归 ，x 为 n 维 列 向 量 ; 


途 为 求 nlinfit 或 nlintool 所 得 的 回归 函数 在 x 处 的 预 
alpha 缺 省 时 为 0.05。 























究 化 学 动力 学 反应 过 程 中 ， 建 立 了 一 个 反应 速度 和 反应 物 含量 的 数学 模型 ， 形 


三 种 反应 物 〈 氧 ，n 戊 烷 ， 异 构成 烷 ) 的 含量 ，/ 








ee 











所 xl an 戊 烷 x2 异 构成 烷 x3 


10 
10 
120 
120 
10 
10 
05 
605 
534 
120 
120 
10 
120 





》 0 


量 ， 并 将 要 拟 合 的 模型 写 








上 nlparci 计 算 回 归 系 数 的 置信 








区 间 ， 月 


昌 nlpredci 计 算 预 测 值 





并 给 出 置信 区 间 。 户 ，… 


























成 国 数 文件 




















$ 2.75 470 
0 14.39 100 
7 2.34 100 
8 4.35 470 
9 13.00 100 
10 8.50 100 
11 0.05 100 


12 11.32 285 
13 3.12 285 
X=X0(:,3:9); 
y=X0(:,2); 





beta=[0.1,0.05,0.02,1.2]; % 回 


[betahat,fj]=nlinfit(x,y，huaxue ,beta); % fj 是 下 面 命令 用 的 信息 


betaci=nlparci(betahat,f]); 











80 10 
190 10 
80 65 
190 65 
300 54 
300 120 
80 120 
300 10 
190 ”120]; 
归 系 数 的 初 值 



































betaa=[betaa,betaci] % 回 归 系 数 及 其 置信 区 间 
[yhat,deltal=nlpredci( huaxue ,x,betahatfj) %y 的 预测 值 及 其 置信 区 间 半 径 ， 置 信 区 间 为 yhat 填 


delta 。 
用 命令 
4 逐步 回归 














逐步 回归 的 命令 是 stepwise， 它 提供 了 一 个 交互 式 画 面 ， 通 过 此 工 


进行 统计 分 析 。 通 常用 法 是 






























































Stepwlse(X，y，inmodel，alpha)， 

















其 中 X 是 自 变量 数据 ，》 是 因 
标 ， 给 出 初始 模型 中 包括 的 子 集 








0.05) 。 


运行 stepwise 命 令 时 产生 三 个 区 形 窗 














Te 












































所 有 这 些 图形 界 面 都 由 热 区 ， 即 当 鼠 标 移 








点 击 后 会 产生 交互 作用 。 








缺 省 时 设 定 为 全 部 


到 


























nlinfit(x,y， Taxue- ,beta) 可 看 到 画面 ， 并 传 出 剩余 标准 差 rmse=0.1933 。 






























































有 具 可 以 








地 选择 变量 ， 

















变量 数据 ， 分 别 为 有 2X 灵 和 7X1 抢 阵 ，inmodel 是 矩阵 的 列 数 指 
变量 ) ，alpha 为 显 间 水平 〈 缺 省 时 为 








: Stepwise Plot，Stepwise Table，AStepwise History 。 

















在 Stepwise Plot 窗 口 ， 显 








示 出 各 项 的 加 














归 系数 及 其 置信 区 间 。 其 中 : 


到 图 形 的 茶 个 区 域 时 ， 鼠 标的 指针 会 变 成 一 个 小 圆 ， 





点 表示 回归 系数 的 值 ， 














点 两 边 的 水 平 〈“ 实 或 虚 ) 直线 段 表 示 其 置信 区 间 《〈 虚 线 表 示 该 变量 的 拟 合 与 0 无 显著 差异 ， 实 
























































线 表 示 有 显著 差异 ) ;绿色 的 线 表示 当前 在 模型 中 的 项 ， 红 色 的 线 表示 当前 不 在 模型 中 的 项 。 
点 击 一 条 线 会 改变 其 状态 ， 即 在 模型 中 的 项 《〈 绿 线 ) 会 被 移 去 〈 变 为 纪 





《红线 ) 会 被 加 入 《〈《 变 为 绿 线 ) 。 次 窗口 


















































[ 线 ) ， 不 在 模型 中 的 项 





PP 的 Export 下 拉 式 菜单 可 以 向 Matlab 工 作 区 传送 各 种 数 


据 。 次 窗口 中 的 Scale Inputs 可 对 输入 数据 的 每 列 进行 正 态 化 处 理 ， 使 其 标准 差 为 1。 








在 Stepwise Table 窗 口中 列 出 了 一 个 统计 表 ， 包 括 回 归 系 数 及 其 置信 区 间 ， 以 及 模型 的 统计 

















量 剩 余 标 准 差 (RMSE)、 相 关系 数 (R-square)、E 值 、 与 FE 值 对 应 的 概率 P。 














例 11 水 泥 凝 国 时 放出 的 热量 y 与 水 泥 中 4 种 化 学 成 分 xLx2x3,x4 有 关 ， 今 测 得 一 组 数据 如 下 ， 
试用 和 逐步 回归 来 确定 一 个 线性 模型 。 




















序号 Xl X2 X3 xx4  y 


Il]2 1 66 9 12 113.3 
]13 10 68 8 12 109.4 
解 编程 如 下 : 
clc,clear 


x0=[L 7 26 6 60 78.5 





2 
3 
4 
4 7 092 0 33 94.9 
0 
孙 
8 
号 


Il10 21 47 4 26 115S.9 
Il1 1 40 23 34 83.8 
Il]2 1 66 9 12 113.3 
]13 10 68 8 12 109.4]; 
X0=(:,2:9); 
yY=X0(:,0); 


StepwWise(X,y) 








得 到 图 Stepwise Plot 和 Stepwise Table 表 。 图 Stepwise Plot 中 四 条 直线 都 是 虚线 ， 说 明 模型 的 显著 





























性 不 好 ， 从 Stepwise Table 表 中 可 以 看 出 变量 x3 和 x4 显 著 最 差 。 移 去 这 两 个 变量 ， 图 Stepwise Plot 
中 点 击 直 线 3 和 直线 4， 这 两 条 直线 变 为 红色 ， 同 时 直线 1 和 直线 2 变 为 实 线 ， 说 明 移 去 变量 x3 和 
x4 后 模型 具有 显著 性 。 从 新 的 统计 结果 可 以 看 出 ， 虽 然 剩余 标准 差 (RMSE) 没有 太 大 的 变化 ， 
但 是 统计 量 F 的 值 明显 增 大 ， 因 此 新 的 回归 模型 更 好 一 些 。 
对 变量 y 和 xl、x2 作 线性 回归 : 
Stepwise(X,y,[1,2]) 
得 到 结果 : 了 =52.5773+1.4683 0+0.6623 2 。 
























































































































































86 建 模 案例 :水塔 水 流量 的 估计 

一 问题 的 提出 

某 居 民 区 供水 机 构 有 一 供 居民 用 水 的 圆柱 形 水 塔 ， 由 于 该 机 构 没 有 测量 流入 或 流出 该 水 塔 
的 水 量 的 装置 ， 水 塔 的 管理 者 只 能 通过 隔 一 段 时 间 测 量 一 次 水 塔 的 水 位 来 估计 水 的 流量 。 但 面 
临 的 困难 是 ， 当 水 塔 水 位 下 降 到 设 定 的 最 低 水 位 时 ， 水 泵 自动 启动 向 水 塔 供水 ， 到 设 定 的 最 高 
水 位 时 停止 供水 ， 这 段 时 间 无 法 测量 水 塔 的 水 位 和 水 泵 的 供水 量 。 通 常 每 天 水 泵 供水 一 到 二 次 ， 
每 次 约 2 个 小 时 。 

水 塔 是 一 个 高 为 12. 2 米 ， 直 径 为 17. 4 米 的 正 圆柱 。 按照 设计 ， 水 塔 水 位 降 至 约 8. 2 米 时 ， 水 
泵 自动 启动 ， 水 位 升 到 约 10. 8 米 是 水 泵 停止 工作 。 
表 1 是 某 一 天 管理 人 员 测 量 该 水 塔 的 水 位 测量 记录 〔〈 符 号 “//” 表 示 水 泵 启动 ) 。 请 您 帮助 
该 管理 人 员 估 计 一 天 中 任何 时 刻 《〈 包 括 水 泵 正在 供水 时 刻 ) 从 水 塔 中 流出 的 水 流量 及 一 天 的 居 
民用 水 总 量 。 




















河 ! 








































































































































































































四 































































































表 1 水 位 测量 记录 

时 0 0.921 1.843 2.949 3.871 4.978 5.900 7.006 7.928 8.967 
水 位 m | 9.677 9.479 9.308 9.125 8.982 8.814 8.686 8.525 8.388 8.220 
时 刻 (D) | 9.9811 10.925 10.954 12.032 12.954 13.875 14.982 15.903 16.826 17.931 
水 位 mm / /10.820 10.500 10.210 9.936 9.653 9.409 9.180 8.921 








时 刻 (b) 19.037 19.959 20.839 22.015 22.938 23.880 24.986 25.908 
水 位 (m 8.662 8.433 8.220 /10.820 10.391 10.354 10.180 














流量 是 单位 时 间 流 出 的 水 的 体积 。 由 于 水 塔 是 正 圆柱 形 ， 截 面 面 积 是 常数 ， 在 水 和 泵 停 上 供 
水 的 时 段 ， 流 量 很 容易 从 水 位 对 时 间 的 变换 率 算出 ， 问 题 是 如 何 估计 水 泵 供水 时 段 的 水 流量 。 
水 泵 供水 时 段 的 流量 与 供水 时 段 前 后 的 流量 有 关 。 如 果 能 够 通过 测量 数据 ， 产 生 若 干 个 时 
刻 的 水 流量 ， 则 计算 水 流量 问题 就 转化 为 插值 问题 。 这 样 水 泵 供水 时 段 的 流量 就 可 以 用 供水 时 
段 前 后 的 流量 插值 得 到 。 作 为 用 于 函数 择 值 的 原始 数据 我 们 希望 水 泵 不 工作 时 段 的 流量 越 准 
越 好 。 
有 了 任何 时 刻 的 流量 ， 就 不 难 计 算 一 天 的 总 用 水 量 。 水 泵 不 工作 时 段 的 用 水 量 可 以 由 测 









































加 | 

















































































































好 



































记录 直接 得 到 ， 如 由 表 1 可 知 从 全 0 到 全 8.967() 水 位 下 降 了 9%.677-8.220=1.457(m)， 乘 以 水 塔 的 
吉 面 积 就 是 这 一 时 段 的 用 水 量 。 这 个 数值 可 以 用 来 检验 插值 的 结果 。 

注 : 此 处 我 们 只 考虑 水 流量 的 插值 问题 ， 当 然 亦 可 以 考虑 水 位 的 插值 问题 ;或 者 亦 可 以 将 
此 问题 按 拟 合 问题 进行 处 理 。 
































































































































三 “模型 假设 

1 流量 只 取决 于 水 位 差 ， 与 水 位 本 身 无 关 。 按照 Torricelli 定 律 从 小 孔 流 出 的 流体 的 流速 〈 流 
量 ) 正比 于 水 面 高 度 的 平方 根 ， 问 题 给 出 水 塔 的 最 低 和 最 高 水 位 分 别 约 是 8.22 米 和 10.82 米 〈 设 
出 水 口 的 水 位 为 零 ) ， 因 为 V10.8278.22 = 工 ， 所 以 可 以 忽略 水 位 对 流速 〈 流 量 ) 的 影响 。 

2 水 塔 中 水 流量 是 时 间 的 连续 光滑 函数 ， 与 水 泵 是 否 工 作 无 关 。 

3 水 泵 是 否 工 作 完 全 取决 于 水 塔 内 的 水 位 的 高 度 ， 且 每 次 加 水 的 工作 时 间 约 为 2 小 时 ， 从 
可 知 两 次 供水 时 间 段 分 别 为 [8.967，10.954] 和 [20.839，22.958]。 

4 水 泵 工作 时 单位 时 间 的 供水 量 大 臻 是 常数 ， 此 常数 大 于 单位 时 间 的 平均 流量 。 

S$ 一 天 之 中 可 以 从 任何 时 刻 开 始 ， 即 从 余 0() 到 人 拓 24() 结 束 与 从 仁 (四 到 关 24+ (结束 差 
别 不 大 。 每 天 同一 时 刻 水 流量 保持 相同 ， 并 且 一 天 中 开始 时 刻 的 不 同 对 一 天 中 的 总 用 水 量 影响 
很 小 。 

6 水 塔 内 部 底面 直径 为 17.4 米 。 
四 “数据 处 理 
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1 水 体积 计算 
广 = 二 万 功 
水 塔 是 一 个 正 圆柱 ， 塔 内 水 的 体积 六 4 ， 其 中 忆 为 水 塔 直径 ， 儿 为 水 面 高 度 。 近 














似 地 取 元 =3.141592654， 得 到 不 同时 刻 水 塔 中 水 的 体积 如 表 2。 














表 2 水 塔 中 水 的 体积 《〈 单 位， 时 刻 《〈 小 时 ) ， 体 积 〈 立 方 米 ) ) 




















时 刻 0 0.921 1.843 2.949 3.871 4.978 5.900 7.006 7.928 8.967 
水 位 2294 2247 2206 2103 2129 2089 2039 2020 1988 1948 

时 刻 9.9811 10.925 10.954 12.032 12.954 13.873 14.982 13.903 10.826 17.931 
水 位 //， /2564 2489 2420 23553 2288 2230 2176 2114 

时 刻 19.037 19.959 20.839 22.015 22.938 23.880 24.986 25.908 

水 位 2053 1999 1948 / 2564 2510 2454 2413 





2 水 流速 (流量 ) 的 近似 
水 流速 度 〈 流 量 ) 是 水 塔 中 水 的 体积 对 时 间 的 导数 。 由 于 没有 水 的 体积 关于 时 间 的 函数 表 
达 式 ， 而 只 是 一 个 离散 的 函数 值 表 2， 因 此 考虑 常用 的 差 商 代 替 导 数 处 理 方法 。 为 提高 精度 ， 采 
二 阶 差 商 公式 。 有 具体 地 ， 表 2 的 数据 点 被 水 泵 两 次 工作 分 割 成 三 组 数据 ， 对 每 组 数据 中 间 数 据 
采用 二 阶 中 心 差 商 ， 两 边 数据 采用 二 阶 向 前 或 向 后 差 商 得 到 数据 表 3。 其 公式 分 别 为 : 
责 2 二 寺 二 庆 3 


12(1 ， 一) 


了 





















































































































































一 Vip+T4a 一 3 


了 

















访 
向 前 差 商 公式 2040H 一 方 ) 
Vs: 3v7 一 人 TV ， 
V 三 
向 后 差 商 公式 2 一直) 





























表 3 水 塔 中 水 的 流速 〈 单 位 : 时 刻 《〈 小 时 ) ， 流 速 〈 立 方 米 /小 时 ) ) 
0 0.921 1.843 2.949 3.871 4.978 5.900 7.006 7.928 8.967 
54.516 42.320 38.085 41.679 33.297 37.817 30.748 38.455 32.122 41.718 


9.9811 10.925 10.954 12.032 12.954 13.873 14.982 13.903 106.820 17.931 














流 // /73.686 76.434 71.686 60.190 68.333 39.217 52.011 56.626 
时 刻 19.037 19.959 20.839 22.015 22.938 23.880 24.986 25.908 
流速 63.023 54.859 353.439 /537.602 57.776 351.891 36.464 








五 ”模型 建立 与 求解 
1 做 出 水 流速 的 散 点 图 


Matlab 程 序 如 下 : 

全 [0 0.921 1.843 2.949 3.871 4.978 5S$.900 7.006 7.928 8.967 10.9$4 12.032 12.934 13.873 14.982 
15.903 10.826 17.931 19.037 19.939 20.839 22.958 23.880 24.986 253.908]; 

I 王 [34.316 42.320 38.083 41.079 33.297 37.817 30.748 38.455 32.122 41.718 73.6086 76.434 71.686 
00.190 608.333 59.217 32.011 30.026 603.023 34.839 33.439 57.602 37.776 351.891 30.464|]; 











plot(tn 十 ); 
title( 流速 散 点 图 ) xlabel( 时 间 (小 时 ) 六 ylabel( 流速 〈 立 方 米 / 小 时 ) ) 
做 出 的 散 点 图 略 。 
2 模型 及 计算 结果 
现在 问题 已 转变 为 根据 流速 VD) 的 一 个 函数 值 表 3， 产 生 函数 y(O) 在 整个 区 间 《〈 二 十 四 小 
时 ) 上 的 函数 或 函数 值 ， 插 值 和 拟 合 是 最 常用 的 两 种 方法 。 如 果 建 立 拟 合 模型 ， 需 要 根据 散 点 
图 的 趋势 ， 选 择 适 当 的 拟 合 函 数 形式 。 若 采用 插值 模型 ， 可 以 考虑 分 段 线性 插值 光滑 性 较 
差 ) 、 分 段 三 次 多 项 式 插值 、 三 次 样 条 插值 等 插值 方法 。 我 们 采用 插值 方法 。 
通过 对 不 同 插值 方法 的 比较 ， 结 合 假设 2， 即 流速 似 时 间 的 连续 光滑 函数 ， 下 面 采 用 三 次 
单条 插值 模型 。 
1) 水 流量 及 一 天 用 水 总 量 ， 首 先 由 三 次 样 条 插值 计算 得 到 水 流量 函数 了 (0D) .其 Matlab 程 序 
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为 : 
估 [0 0.921 1.843 2.949 3.871 4.978 5.900 7.006 7.928 8.967 10.954 12.032 12.954 13.875 14.982 
13.903 16.826 17.931 19.037 19.959 20.839 22.958 23.880 24.986 25.908]; 
I=[34.316 42.320 38.083 41.079 33.297 37.817 30.748 38.433 32.122 41.718 73.086 76.434 71.080 
00.190 68.333 59.217 52.011 56.026 63.023 354.859 35.439 357.602 57.776 51.891 36.464j]; 


[1Ln]=size(b; dl1=t(n)-t(1); 








x=t(0):1/3600:ttm); % 被 插值 点 














ys=interpl(bPx，spline ); % 样 条 搬 值 输出 


plot(Cxys); 





titleC 样 条 插值 下 的 水 流速 图 ); xlabel( 时 间 〈 小 时 ) )iylabel( 流速 〈 立 方 米 /小 时 ) ) 
做 出 的 水 流速 曲线 图 略 。 
2) 一 天 (24 小 时 ) 的 用 水 总 量 ， 用 三 次 样 条 插值 模型 得 到 的 函数 /iD) 在 区 间 [0.24] 上 积分 
得 到 的 结果 为 1257.3 立 方 米 。 
3) 误差 分 析 直接 由 表 3 的 水 流速 数据 ， 用 梯形 公式 进行 数值 积分 得 到 的 结果 为 12$0.3 立 
方 米 ， 与 上 面 得 到 的 结果 比较 ， 则 有 : 三 次 样 条 插值 的 绝对 误差 为 7 个 立方 米 ， 相 对 误差 为 
0.56% 〈 不 到 百 分 之 一 ) 。 
另外 ， 测 量 数据 的 误差 可 控制 在 0.5%; 数值 微分 、 数 值 积 分 的 误差 亦 可 以 控制 在 一 定 的 范 
围 内 。 
六 ”模型 的 稳定 性 分 析 与 检验 





































































































































































































用 不 同时 刻 作为 起 始点 ， 使 用 三 次 样 条 插值 函数 得 到 的 水 流量 函数 /0) ， 在 长 度 为 24 小 时 
的 时 间 区 间 上 积分 ， 所 得 结果 相差 无 几 ， 详 见 表 4。 这 说 明 我 们 所 建立 的 数值 微分 模型 和 三 次 样 
条 插值 模型 是 稳定 的 。 





















































表 4 不 同 起 点 计算 出 的 24 小 时 用 水 量 〈 单 位 : 立方 米 ) 














起 始点 0.00 0.25 0.30 0.753 1.00 1.253 1.30 1.75 








2 检验 

1) 分 三 段 (水泵 未 工作 ) 的 实际 用 水 量 与 由 模型 推算 的 用 水 量 的 差异 很 小 ， 见 表 5， 最 
大 的 不 超过 4%， 三 段 总 计 不 超过 0.5%。 其 中 由 模型 推算 的 用 水 量 为 水 流量 函数 /2) 的 积分 ; 
实际 用 水 量 为 水 塔 内 水 的 体积 之 差 。 




















































































































表 5 分 三 段 的 实际 用 水 量 与 模型 计算 出 的 24 小 时 用 水 量 比较 〈 单 位 : 立方 米 ) 
国 和 | 量 ”模型 用 水 量 “绝对 误差 “相对 误差 
第 一 段 [0，8.967] 345.3682 “338.0$11 7.3171 ”2.1% 























第 二 段 [10.954，20.839 616.3161 ”620.8341 ”4.5180 0.7% 
第 三 段 [22.958，25.908 151.7308 ”156.6000 ”4.8692 ”3.2% 
1113.4151 ”1115.4852 ”2.0701 ”0.2% 














2) 两 次 充 水 期 间 ， 水 泵 注入 量 的 差异 大 约 为 3 个 立方 米 ， 不 到 0.5%。 
水 泵 充 水 量 = 充 水 后 的 水 量 + 充 水 期 间 的 流出 量 一 充 水 前 的 水 量 。 
第 一 次 水 泵 充 水 量 =2564+117.4463-1948=733.4463 立 方 米 ; 
第 二 次 水 泵 充 水 量 =2564+120.7052-1948=736.7052 立 方 米 。 
于 水 和 泵 功率 = 充 水 量 / 充 水 时 间 ， 由 上 述 计算 知 ， 水 泵 的 功率 大 约 为 367 立 方 米 /小 时 ， 而 且 两 
次 冲 水 期 间 计 算出 的 功率 大 致 相等 ， 这 与 实际 问题 是 一 致 的 。 
3) 用 水 高 峰 的 比较 

















































































































实际 与 模型 之 间 相 差 无 几 。 实 际 用 水 高 峰 可 近似 地 用 差 商 最 大 值 点 表示 为 交 11， 即 上 午 11 
点 钟 左 右 ; 模型 得 到 的 用 水 高 峰 可 由 模型 得 到 的 水 流量 函数 依次 在 单位 区 间 上 积分 或 者 找 出 水 
流量 函数 的 最 大 值 点 得 出 ， 它 们 都 在 11 点 左右 。 
七 ”模型 的 优 缺 点 及 推广 
1 优点 
(1) 模型 灵活 性 好 、 稳 定性 强 ， 可 用 于 那些 拥有 地 方 性 的 竖 直 圆柱 型 水 塔 的 小 城镇 和 乡镇 。 
模型 中 的 输入 数据 可 以 是 任何 近似 均匀 的 时 间 间 隔 时 的 水 位 ， 时 间 间 隔 大 约 2 小 时 。 
(2) 模型 中 的 数学 概念 简单 ， 并 且 容 易 理 解 。 只 用 到 数值 计算 知识 。 
43) 模型 容易 实现 ， 且 给 出 了 一 天 里 水 流速 度 和 用 水 量 的 精确 估计 。 
2 缺点 
(1) 本 模型 受 水 塔 的 几何 形状 限制 。 
《2) 光滑 曲线 的 逼近 方法 不 能 模拟 真实 流动 中 流速 的 微小 变化 ， 实 际 流 动 中 流速 可 能 会 
种 高 程度 的 “噪音 ”， 即 激 波 出 现 。 
3 改进 与 推 
(1) 我 们 可 以 在 模型 中 用 一 个 参数 来 限定 不 同 几何 形状 的 水 塔 。 
《2) 可 以 通过 对 流速 数据 进行 回归 分 析 等 一 系列 处 理 ， 以 便 得 到 一 些 随 机 变化 的 特征 。 





























































































































































































































































































































， 习 题 十 
1， 用 给 定 的 多 项 式 ， 如 了 =X 一 6x 十 5 一 3， 产 生 一 组 数据 (0， 芒 ) ， 守 = 12 ， 
再 在 世上 添加 随机 干扰 〈 可 用 rand 产 生 〈0, 1) 均匀 分 布 随机 数 ， 或 用 randn 产 生 人 (0 分 布 
随机 数 ) ， 然 后 用 于 和 添加 了 随机 干扰 的 六 作 3 次 多 项 式 拟 合 ， 与 原 系数 比较 。 如 果 作 2 或 4 次 
多 项 式 拟 合 ， 结 果 如 何 ? 
2， 已 知 平面 区 域 0&xY< 3600 ，0 和 和 4800 的 高 程 数据 见 下 表 〈 单 位 : m) 。 


1410 960 940 880 800 430 210 
1440 1140 1110 1050 950 540 300 
1470 280 1080 620 370 
1500 510 1300 850 550 
1550 1 1600 1500 1550 




















1200 1550 900 1150 





1200 1 1200 880 1050 








1500 1100 870 936 
1450 1 700 810 780 
1250 780 650 
1020 500 480 
850 200 350 
670 300 150 
了 /和 友 10 400 800 1200 1600 2000 2400 2800 3200 3600 4000 4400 4800 5200 5600 


试用 二 维 插值 求 汪 了 方向 间隔 都 为 50 的 高 程 ， 并 画 出 该 区 域 的 等 高 线 。 
















































































二 Dx _ 中 
3， 用 最 小 二 乘法 求 一 形 如 上 二 “5 的 经 验 公 式 拟 合 下 列 数据 
1 双 3 4 5 6 7 8 
15.3 20.5 27.4 36.6 49.1 65.6 87.87 117.6 




















4.， 某 人 记录 了 21 天 使 用 空调 器 的 时 间 和 使 用 烘 干 器 的 次 数 ， 并 监视 电表 以 计算 出 每 天 的 
耗 电 量 ， 数 据 见 下 表 ， 试 研究 耗 电量 (KWH) 与 空调 器 使 用 的 小 时 数 (AC) 和 人 烘 干 器 使 用 次 
数 (DRYER) 之 间 的 关系 ， 建 立 并 检验 回归 模型 ， 诊 断 是 否 有 异常 点 。 

序号 123453 6789 10 1 

KWH 35 63 66 17 9%4 79 9%3 66 94 82 78 













































































AcC 1.3 4.3 3.0 2.0 8.3 6.0 13.3 8.0 12.3 7.3 6.5 
DRYER 1 2 203 31 1 1 2 3 

序号 l]2 13 14 13 16 17 18 19 20 21 
KEWH 663 77 73 602 85 43 37 33 .63 33 

AC 8.0 7.3 8.0 7.3 12.0 0.0 2.3 3.0 7.3 0.0 
DRYER 1 2 2 1 10 3 0 1 0 




















5， 在 一 丘陵 地 带 测量 高 程 ，x 和 y 方 向 每 隔 100 米 测 一 个 点 ， 得 高 程 如 下 表 ， 试 拟 合 一 曲面 ， 
确定 合适 的 模型 ， 并 由 此 找 出 最 高 点 和 该 点 的 高 程 ， 给 出 回归 分 析 。 
X 












































了 100 200 300 400 500 

100 636 697 624 478 450 

200 698 712 630 478 420 

300 680 674 598 412 400 

400 662 626 552 334 310 

6. 一 矿脉 有 13 个 相 邻 样本 点 ， 人 为 地 设 定 一 原点 ， 现 测 的 各 样本 点 对 原点 的 距离 x， 与 该 

单 本 点 处 某 种 金属 含量 y 的 一 组 数据 如 下 ， 画 出 散 点 图 观测 二 者 之 间 的 关系 ， 试 建立 合适 的 回 
归 模型 ， 如 二 次 曲线 、 双 曲线 、 对 数 曲 线 等 。 

x23 4 5357 8 10 

y_106.42 109.20 109.58 109.50 110.00 109.93 110.49 

















下 
























































X ll 14 113 16 13 19 21 
110.39 110.60 110.90 110.76 110.00 111.20 110.33 
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